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Preface 


This material is designed for students Interested In the area of space flight 
planning, who after training, may serve as flight planning aides. Flight plan- 
ning aides are Individuals who will assist engineers and scientists In the plan- 
ning c ; orbital space missions. With the development of the NASA space shut- 
tle ar>u its earlier successful missions, the time will come shortly when orbital 
flight missions will become a routine and nearly continuous activity. Space 
station construction Is likely to be the first significant activity of the shuttle, 
along with the deployment of strategic and scientific experiments. Extensive 
commercial use of orbital flight is expected to follow. 

With orbital flight becoming much more common, the need for routine flight 
planning activities increases. There are significant tasks In the planning of 
flights that can be performed by individuals with less than ari engineer's 
background. Those tasks will be performed by the flight planning aide. 

The objective In the training procedure for a flight planning aide Is to provide 
an output, In a short period of time, that is capable of performing engineering- 
type calculations and analyses. To perform useful and significant calcula- > 
tions, it is not necessary that one know the details and background of all 
operations. That Is true for tasks other than calculations. While it Is hoped 
that calculators and computers do not retard us in ability to add, multiply and 
reason, they do allow many of us who bore to error with the task of addition to 
become painlessly highly accurate. Some of us are curious about process and 
insist upon a thorough understanding. The flight planning aide student will be 
capable of performing certain operations without understanding their bases, 
and hopefully have a curiosity that will be required for learning extension, k is 
on that basis that this training material has been put together in a directed 
ler-nlng approach. 
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For this material, a calculator Is necessary. Calculators and computers, 
along with the pencil, will certainly be the professional tools of the flight plan- 
ning aide. Use of the calculator, common now In the genera! population. Is em- 
phasized first In the workbook In order to proceed quickly toward programming 
and eventual use of the computer. A strong objection might be made on the 
emphasis put In handling a particular calculator, the Hewlett-Packard 67. That 
we have done. Any progiammable calculator with a good guide to Its use, 
however, will be appropriate for the given exercises. Texas Instruments has 
such a calculator, though It uses algebraic logic, and there are others. 
Although we have described the HP keyboard and made specific reference to 
the HP handbook, any calculator may be used with the understanding that Its 
operations must be thoroughly un<ierstood, through use of Its calculator guide 
or owner's handbook. The first requirement for the flight planning aide stu- 
dent, therefore. Is to be able to learn the use of the calculator. A strong 
mathematical Inclination Is also thought necessary. 

But we have attempted to Integrate learning activities. Flight planning 
vocabulary and orbital flight relations are Introduced In the beginning and are 
used while facility with problem solving and the first major objective, use of the 
calculator. Is learned. To do this as we have, primary calculator instruction 
must come from the manufacturer or elsewhere. The learning tasks proceed 
through the solving of problems, as relatively more complex problems are 
quickly Introduced Into the calculator and programming exercises. 

Teaching beginning students through such a problem sequence has long 
been a fascination of mine, but difficult to effect In traditional college sequen- 
tial course curriculum. An Information scattering technique Is used, one I use 
In lecture and believe to be quite useful In Inputing data to the learner that Is 
not to be Immediately used. But too much non-useful written Information can 
be Inhibitive. This written material contains only purposeful information with 
the hope that curiosities will be aroused to the point of asking questions and 
seeking further knowledge. The teaching-learning philosophy of Mr. Davis, in 
teaching such a course in flight planning, has been quite similar to my own. 

H.E. Blackwell, Ph.D. 

June 30, 1981 
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A. HAND CALCUUTOR USE AND 
INTRODUCTION TO PROBLEM ANALYSIS 


A.O. INTRODUCTION 

The flight planning aide is expected to solve many problems inwolvinq 
numerical calculations. Some problems m. jht be done easily using paper and 
pencil algorithms and simple baste computations. Other problems will be 
mure difficult to i»nalyze and calculate, and may be handled most efficiently 
through use of a calculator. Since high accuracy will be necessary, calculator 
use allows lengthy, cumbersome calculations to be handled with speed and 
accuracy. Take full advantage of the power of the calculator. 


A.1. CALCUUTOR LOGIC AND SIMPLE OPERATIONS 


Three types of logic are used on calculate, s; Reverse Polish Notation (RPN), 
Algebraic Logic (alg) and Arithmetic Logic (Arith ). Thos e calculators using 
RPN are characterized by a key marked flinter |or [~Sive^ . The algebraic and 


arithmetic calculators are distinguislred by the order in which some operations 
are done. The key sequences below show how each type of calculator would 
be used in a simple operation such as: (4 + 5) x 3. 


RPN Alg Arith 

4 4 4 

enter f + + 

5 5 5 

+ = X 

3x3 

3 s 

= I ANS 

ANS 

These notes are written for a Programmable Calculator which uses HPN. 


X 


ANS 27 


1 





The Hewlett Packard HP-67 has been chosen, Note two conventions we will 
follow; operations will be written in calculator language, not In algebraic 
language, with the language written in a column, top to bottom. These; conven- 
tions i^re motivated by practical considerations; the calculator language gives 
preparation for later programming activities, and the column representation 
enables you to write mnemonic notes on the same line to the right of each 
operation. The RPN calculator was chosen in order to de-emphaslze the 
algebraic manipulations and to emphasize the calculator language. 

Your initial task will be learning the use of various functions and 
operators on the calculator. Read the material from pages 15-24 in the 
HP67 Owners Handbook and Programming Guide to get a “ieel" for the 
calculator. Do the suggested exercises. These will provide some sense of 
the power of the calcijlator In performing simple calculations as well as 
some Idea of what is involved In writing and running a program. Read the 
materials from Part 1, section 1, pp. 27-M. Note the Information concern- 
ing the keyboard on page 27. Note further that each key on the keyboard 
can perform four different functions. One function Is indicated on the flat 
plane of the key -e; another Is printed In black on the slanted (ace of the 
key; a third and fourth function is indicated by printed symbols in gold and 
blue below the key. There are three prefix keys; a gold one, a blue one, 
and a black one. By pressing one of the prefix keys before pressing a 
function key, you select the function printed on the slanted key face or 
one of the functions in gold or blue below the function key. As you read 
the material In section one, work through the exercises. If the need for ad- 
ditional problems is felt, do the following; 

(14 - 3) + [(8 X 5) 4] ANS; 21 

[ -2 -r 4 (1- 6)1 _2 

11 

Y(4 - 2) -!■ (15 + Sj -r ■V(6 - 5) -f- (4 X 20) ANS; 14 

In evaluating an expression with pencil and paper involving several opera- 
tions, 'the following procedure is conventior\ally followed; 

(1) Perform the operations within each grouping symbol, if 
any, beginning with the innermost grouping symbol and 
working out to the entire expression. 

(2) Perform multiplications and divisions in order from left 
to right. 

(3) Perform additions and subtractions in order from left to 
right. 
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In using tne calculator, a decision rnust be made on where to begin before 
the first Key Is pressed, The above procedure will help you. Please note 
that in general, only one operation Is performed at a time, 

Now proceed to read and work through section 2 on Display Control (pp. 
41->51), This section explains the methods for displaying numbers In certain for- 
mats, Including scientific notation, fixed point display, and engineering nota- 
tion display, Work through these exercises also. Read f' e Information on 
automatic switching between fixed and scientific notation If the answer Is too 
large. Do the related exercises. Note the information in this section relating to 
"Keying In Exponents of Ten,” "Calculator Overflow,” "Error Display,” and 
“Low Power Display,” Remembering where certain types of Information may 
be found in the handbook will facilitate Its use as a quick reference. 


A.^ USE OF THE ,;\UTOMATIC MEMORY STACK 


Th'S section of the HP-67 handbook deals with the four register 
automatic memory stack and the LAST X register of the HP-67 Calculator. 
Thero are four registers In the calculator that are positioned to form the 
automatic memory stack. These registers are labeled X, Y, Z, and T, They 
are "stacked” one on top of the other with the displayed X-regIster on the 
bottom, When the calculator is switched on, these four registers are 
cleared to zero, and the display shows 0.00 In RUN mode. This represents 
the contents of the X-reglster. Certain keys on the calculator enable the 
stack contents to be reviewed or to shift data within the stack for com- 
putation at any time, 


When the h function key (black) and the ( Ri ) key (roll down) are pressed, 
the s tack contents shift down one regi-ster. If the h function key and the 
1 1 |(roll up) key are pressed, the stack contents roll up instead of down. 

r"x*Tvi 

Pressing the h function key and the ^ ' (x exchange y) key causes the 


contents of the x and y-registers to interchange without affecting the Z 
and T registers. Pressing the h function key and the 


key a second 


time In succession will restore the numbers in the X and Y register to their 
original position. 

The g function key (blue) and the STK | key will allow the entire stack 
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contents to be reviewed at any time, This is particulariy useful in checHino 
to be sure that the pr oper numbers have Iwen enteied into the register, 
When operating the g 


STK 


key manually, y ou can slow down or speed up 

the review of the stack contents by pressing or any other key on the 
keyboard while the calculator Is executing a g F^tKl review. 

The fcO n (clear X, key) is used to clear the x-register to zero. No other 
register Is affe cted when you press clear X, When you are entering several 
numbers, the f ENTER tl key on the calculator Is used to separate the 


digits of one number from another, 

In addition to the four strck registers, there Is a last X [Tsf x1 register, 
used wltit the h function key. This key Is used to presen/e the value that 
was In t he displayed X register before the performance of a function. The 
1ST X key makes it easy to recover from keyboard mistakes such as 


pressing? the wrong function key or keying In the wrong number. The last X 
LST X key is also useful In calculations where a number occurs more 
than once. 


In the past few paragraphs, we have reviewed very briefly material covered In 
pages of the Owners Handbook. Read these pages with calculator in 
hand and do the suggested exercises while reading. 


A.3. EVALUATING FORMULAS 

A fornula is a statement In mathematical language of a general rule or 
principle. It tells how the value of one unknown may be found by perform- 
ing certain operations upon the variables whose values are known. The 
unknown wilt generally be placed to the left of the equality sign such that 
it fs evaluated by performing the operations on the right. 

unknown = [variable expression] 

A calculation to determine a numerical value for the unknown thus in- 
volves (a) replacing the variables with their numerical values, and (b) per- 
forming the indicated operations of addition, division, etc,, to arrive at 

unknown = number 

In the calculator exercises which follow, fomiulae used In astro- and aero- 
dynamic calculations will be presented in that format, unknown on the 
left. Variables will be defined. While a clear understanding of the formula 
and Its basis is not necessary at this time, some understanding of the 
many terms and relationships will follow through a careful reading and 
performance of the exercises. 
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In addition, constants used In the actual plannlr,g of missions will he 
used for the exercises In order to facilitate use of accurate and consistent 
numerical valuee. For your convc;'.!sr,ce, data from the Astrodynamical 
Constants Handbook, JSC-14^, is t«pii)duced as Appendix A. The con- 
stants are written In powers of 10 notation, For example, the value of the 
ofsrth’s equatorial radius Is presented as 


1 E.f. = .2092572178477090 0(» ± .16404 002 Int. K. 

The -^006 means that the l&digit number is multiplied by 10 to the 8th 
power (10*). The ± .16404 -f 002 provides the accuracy of the given 
number. Written without power of 10 notation, the above constant would 
be 


1 E.r, = 20, 925,721.78477690 ± 16.404 int. ft. 

Plus or minus 16 ft in a measure of 21 million ft is not bad accuracy, 
Strictly speaking, because 16,404 is the possible error, It Is impractical to 
maintain more than three decimal places here in the measurement. Note 
that the units of measure are given and check the Appendix to see that the 
radius is given in several other units of measure. It Is important that ap- 
propriate an<'* consistent units are used. We will discuss units and their 
conversion IfAner, To Input powers of 10 to the calculator see pages 48*51 
of the Owneis Handbook. 

As a first example, consider the following information. The acceleration 
of gravity at the earth’s surface is given by the formula 


where 

fji Is the gravitational constant, In 
units of international feet per 
second (Int. ft/sec’) .squared 

^ Is the earth's equatorial radius, In int. ft. 

Taking g as the unknown, your exercise is to calculate its numerical value. 
The calculator steps or program sequence is as follows. 
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0.200257721786 E 08 


Int. H. 


t 




gx 


[~a}^ 

These are the cal'^ulator steps for evaluation of the variable expression on 
the right side In the given formula. To carry It through, numerical values of 
n and Rflp must be obtained. They may be obtained, In appropriate units, 
from Appendix A. While you may enter the numbers directly into the 
calculator. It will be helpful to first list those numbers with their units to 
the right of the variable In the given sequence. We have provided B,,. 
Compare your answer with that of g as jjlven In Appendix A. 


A.4 USE OF ASTRODYNAMiCAL CONSTANTS 

For evaluation of any given formula, dimensions or units used in the for- 
mula must be consistent. If two of the variables are length, for example, 
one cannot use nautical miles (n. mi) for one and international feet for the 
other. A decision must be made as to which of the units is to be used. If 
we have a numerical value in one set of units only, it may be converted to 
other units by use of a conversion factor. While conversions will be 
discussed later, it is generally simpler to use a multiplicative factor. For 
example, multiplying n. mi by a constant factor equal to 6076.1155 con- 
verts the number to Int.ft. To convert ’nt.ft. to n.ml we multiply by the 
reciprocal of that number. 

Exercise A.4.1. Convert the earth's equatorial radius (R.^) in 
nautical miles (n.mi) to the radius in international feet (int.ft). 

R«, n. mi 

t 

Conv 6076.1155 IqL«- 

n.mi 
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Int, ft, 


ANS 


A given exercise or problem consists of a set of words from which you 
must extract the set of "given" Information (Input) which relates to for- 
mulae to be used, and determine what you are asked to "find*' or calculate 
(output), Relate quickly the given Information to formula variables and 
check the units for consistency. Calculator steps are provided In the In- 
itial exercises. 


The altitude of a satellite Is the height above the earth's equatorial 
radius. For a satellite In a circular orbit, the altitude Is the same 
everywhere in the orbit, by definition (see Figure 1A), The actual height 
varies since the earth is not truly a sphere, but we will discuss this later, 
The formula for finding the semimajor axis (equivalent to radius here) for a 
satellite in circular orbit Is 

a = H (circular orbit) 

where 

a Is the semimajor axis, in units Identical to units of H and R^ 

H is the altitude, In units identical to R,«, 

Is the earth’s radius, in units of nautical miles (n.ml) 

Exercise A.4^ If the shuttle orbiter is in a circular orbit of 160 n.ml 
altitude, find the semimajor axis. 

H 160 n.ml 

t 

R«, n.ml 


I ANS I 

Find the semimajor axis In Int.ft. 

The velocity of a satellite in a circular orbit is everywhere the same, and 
is given by 


c. 
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Vok« = 


where 

a is the semimajor axis, 

M Is the gravitational constant. 

Exercise A.4.3. Find the velacity of the orbiter In ft./sec. when it is In a 
circular orbit of 160 n. ml. 


H 

160 

n, ml 

t 





n. mi 

Conv 

.60761155 

lnt.fi 

n.rfll 

EEX 




4 

X 

H Int. ft.3 

sec“ 

g xs y 
ANS 1 


Note that at one point In the above calculations, the semimajor 
axis In International feet Is displayed. Does this value agree with 
your answer in the previous example? 


The time required to travel between points Is given by 

t = 1. 

V 
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where 

t is the time lapse, In secs, 
d Is the distance In Int.ft 
V is the velocity In Int.ft per sec. 

Exercise A.4.4. Tf>e distance between a tracking site near Los 
Angeles and a point near Baytown, Texas Is 1210.83 nautical miles. 
How long will It take the orbiter to traverse an equivalent distance 
If It is In a circular orbit of 160 n.mi altitude? (see exercise A.4.2.) 


d 

1210.83 

n.ml 

t 



Conv 

0.60761155 

Int.ft 

n.mi 

EEX 




4 

X 

V from previous exercise / Int.ft 

V sec 


ANS 


Using the calculator function, express this answer In hours, 
minutes, and seconds. 

The total distance around a circular orbit is equivalent to the cir- 
cumference of a circle of radius a. 

d = 2 rr a 

where 

d is the distance 
a Is the semimajor axis 

Exercise A.4.5. If the orbiter is in a circular orbit of 160 n. mi. altitude, 
how long does it take to complete one revolution In this orbit? 
(Remember that T -= d/v.) 
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2 


I 

g n 

X 

a 

X 

V 

+ 


ANS 


A.S. THE “EXPONENTIATION”OPERATOR 

The operation fy" has many uses; it can be used to find squares, cubes, 
square roots, cube roots, and so on. It Is also used to raise a number in the y 
register to a fractional power. 

The period of a satellite in an earth orbit is the time required to make 
one complete revolution in the orbit. The period is given by the formula 
(the orbit need not be circular) 


where 



P is the period. In seconds 
M is the gravitational constant, in consistent units 
a is the semimajor axis, in consistent units. 


Exercise A.5.1. What is the Period of the orbiter when in a circular orbit 
of 160 n.ml. altitude? 

a 160 n.mi 

f 

Rgq n.mi 
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+ 

3 

f 

2 

+ 

f/ 

gjT 

X 

2 

X 


^ CL-iSJ! 

sec® 

t i~~% 


ANS 


Compare this answer with the one you found in exercise A.4.4. 


Exercise A.5^. Suppose a communications satellite is in a circular 
orbit with an altitude of 19322.91 n. mi. What is the period in hours, 
minutes, secs.? 


H 19322.91 n.ml 

t 

Roq n.mi 


1.5 

f/ 
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grr 


X 

2 

X 

M 

iVx 




Conv 


n.ml’ 

sec* 


3600 


sec./hr. 


f fix4 
f— H.MS 
~ANfl 

If this satellite is directly over the equator south of Houston at mld- 
nite, v\/here will it be approximately one hour later? 


The time required for a rotating body to rotate through a given angle is 
given by the formula 



where 

t Is the time lapsed, in seconds 

d is the total angle traversed, in degrees or radians 

u) is the angular velocity, in degrees per sec or radians per 
sec. 


Exercise A.S.3. The angular velocity of the earth relative to a 
processing equinox given in Appendix A is 

^ p = .729211585 E - 04 radians per sec 


13 


Use the Appendix to find this vaiue, Find the time, in hours, 
minutes, seconds, required for the earth to rotate 2 it radians 
(equivalent to 360*). 

on 

t 

2 

X 

p radfsec 

Conv — 1 — hr/sec 

3600 


X 

f FIX 4 
f — H.MS 


ANS 


Compare this ans\wer with the period of the communications 
satellite. 


A.6. “READING” EQUATIONS AND PROBLEMS 

We will have a distinct advantage when we are able to "read" an equation 
and translate it Into caiculator language with the same ease that we read 
alphabetic symbols and translate those into English. Some simple rules are 
helpful: try to do terms within multiple parentheses “Inside out” taking advan- 
tage of the temporary storage features of the shift registers. In very long com- 
plex expressions, it may be neccessary to store intermediate results in one of 
the working storage registers provided by the calculator. 


Exercise A.6.1. Practice order of computation using the following. 
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X = (14* + 196'^* ) (Vi) 

X + (72/9)'^’ - 5 

x= 1+3 (6 4 2(7 + 13(7))) 


But a major objective for us is to be able to solve problems related to flight 
planning. So definitions, concepts, principles, and ideas presented here will 
relate directly to actual problems that are encountered in flight planning. To 
assist In problem "reading,” consider four phases that are Involved In problem 
solving, First is an understanding of the verbal statement of the problem. 
Understanding the problem Involves being able to answer these questions; 
What data is given? What are the unknowns? What conditions are given? For 
many this can be difficult. 

The second phase In problem solving Involves seeing how the various Items 
are connected, how the unknown Is linked to the given Information In order to 
obtain the idea of a solution. This involves making a plan. You have a plan 
when you are able to determine the appropriate equations or relations and 
outline the calculations, computations, or constructions that will have to be 
made In order to arrive at a solution. Here, your plan will Invol.n primarily the 
Identification of the appropriate formula(s) relating the unknowns lO the given 
data and determining the calculator process (step-by-step) necessary for 
evaluation of the unknown through use of the formula. Developing the plan is 
the main part in the solution of the problem. 

The third phase involves carrying out the plan, that Is to say, performing the 
calculations. Since the calculator will do the work, this is the simplest part. 
The last step Is to look at the completed solution and see if it Is reasonable. 
Does It satisfy the conditions of the problem? 

Consider now the set of formula given below. They are used In determining 
performance characteristics of a rocket engine. 


(1) 

'•P = F/Wb 

(11) 

F = Wbl^ 

(III) 

Wb = F/l^ 

(Iv) 

U - c/g 

(V) 

c = gi^ 
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where 


F is the rocket thrust In pounds (lbs); 

Wb is the fuel burned per ^ec, Ibs/sec; 

c is the exhaust velocity of the gases issuing from the rocket 
nozzle, in ft./sec; 

l,p is the thrust developed for each lb . of fuel burned In one se- 
cond, 1^ Is called "the specific lmi?ulse” and is a property of 
the fuel used; 

g Is the acceleration of gravity, In ft/sec’, 


Exercise A.6.2. Each S R B at orbiter launch burns fuel at the rate of 
9,330 Ibs/sec. The l,p of the rocket fuel used Is 260 secs. 

a) What is the thrust of an SRB? 

b) What is the exhaust velocity of the burning gases? 

The nain engine (SSME’s) each burn a Hydrogen-LOX fuel at the rate 
of 1 <320 lbs. per second. Each develops a thrust of 470,000 lbs. 

c) What is the l,^, of the Hydrogen-LOX fuel? 

d) What is the effective exhaust velocity of the burning 
Hydrogen-LOX fuel? 

Give the dimensions for each answer. There are 2 SRB's and 3 SSME’s. 
The total weight of the SRB's, orbiter, and external tank at launch Is 
4.461485 E06 lbs. 

e) What is the total thrust at launch? 

0 Will the system get off the pad? 


Let us review the exercise just completed. The first two statements provided 
you the following: 

Wb = 9,330 Ib/sec 
l,p = 260 sec 

You related the words "burns fuel at the rate of” with the definition of W^ given. 
To answer the first question, you looked at the symbol for thrust given in the 
definitions and determined immediately that formula (II) was to be used. From 
the formula 

F = Wp l,p 
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a slmpie multiplication determines the answer. 


Ib/sec 


Wb 

t 

X 


9,300 


260 

2,418,000 


sec 


lb 


For question (b), formula (v) was determined appropriate and cejculation pro- 
ceeded similarly, To answer question (0> it was necessary to conclude that the 
total thrust must be greater than the total weight for lift off, a reasonable con- 
clusion one might reach without an understanding of Newton’s laws of 
physics. 

Note that In the set of five formula, there are but two primary relations. The 
other three are simply equivalent forms that we have provided you In order that 
the unknown variable be positioned on the left. For example, the first was 

I =: L 


Multiplying both sides of the equation by W), yields 


IspWb - ^ 


Cancelling Wp on the right side isolates F: 

'.pWb = F 

This is the same equation, given In what is called an equivalent form. Now if 
we divide both sides of the equation by l^,. 


XWb = _F 

i»p 


or 


Wb = 


F 


The variable is isolated giving another equivalent form of the equation. We 
use this relation when Wt, is the unknown variable to be calculated. To simplify 
your work in many cases, we will provide both the basic equation or formula 
and the useful equivalent forms. 
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A.7. DIMENSIONS AND THEIR CONVER^N 

Numbers obtained from most physicol measurements are of two types: pure 
numf^rs, and numbers with dimensions of mass, weight, and time, Pure 
num'bers are nearly aiways formed as ratios of dimensional numbers, so that 
dimtnsions cancel out in the resuit. Mechanical measurements yield numbers 
with named dimensions or units that are combinations of the basic dimen- 
sions or units or mass, iength, and time. Tor example, the unit for force is the 
pound, but the basic are mass X iength/time^. When evaluating a formula, we 
should perform a dimensional analysis If the final units are in question. For ex- 
ample, in the last exercise, dimensional analysis of the formula F = W„l,p 
would yield 

X 1^ (ped) = F (lbs) 

A conversion constant has a riiitio of dimensions associated with it. Obtain- 
ing a conversion constant is the first step In converting dimensions. Let us 
take time as an example. There are, of course, si:dy seconds in one minute, or 
we nay say that 60 seconds is equivalent to 1 minute. 

60 secs - 1 min. 


Dividing both sides by 1 min, we obtain 


60 secs _ 1 j»kT 
1 min ~ 1 min 


60 secs 
min 


Note that on the -Ight is the pure number 1, and on the left we have obtained a 
conversion constant. Since multiplication by 1 does not change a number, we 
may use the conversion constant as a multiplicative factor. Note also that had 
we divldel by 60 secs, we would have obtained 

•j _ 1 min 
60 sec 

the reciprocal of the above factor, To convert seconds to minutes, j multiply 
using this factor. To convert minutes to seconds, we use the other. 
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As an example, we are asked to convert to minutes t * 14.86 seconds. The 
object is to use the conversion constant with seconds In the denominator so 
that it will cancel the seconds in time t, 

t = 14,85^ X -Kia = 0.270 min 
60 sec 

Or, suppose we ana to convert to ft/min the velocity v 14.85 ft/sec. We do the 
following; 


V = 14,86 X 60 ^ = 891.00 -4- 
.sac min min 

We have used a conversion constant in previous exercises, to convert nautical 
miles to international feet. From the data fin the section Equatorial Earth 
Radius) in Appendix A, we had obtained 

1 E.R. = .2092572178477690 + 008 Int.ft. = .344393088K19158 + 004 n.ml 

tsSaoSsSigiME S 'nj^' ^ ^ 

.60761 1549E 04 - 1 

n.m( 

Or, there are 6076.1155 Internationa! /eet in one nautical mile. Check also In 
Appendix A the section on "Equivalents and Conversion Factors." 

Control of dimensions will be made easier If we multiply when converting 
and check to make certain that the dimension to be converted from always 
cancels out in the multiplication. 

• 

The earth's rotational velocity is 

W p = ,72921 1585E - 04 rad/sec 

There are 2 n radians in one revolution. There are 86,400 seconds in one day. 

Exercise A.7.1. Find the angular velocity of the earth in revolutions per 
day. 

0.7292585E - 04 Ta^ X ± rev x 86,400 Seg 
'9eg 2 rr Tad. 1 day 

= 1.0027379 rev/day 
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How many days per revolution? If there are 24 clock hours per day, 
how many hours per revolution? 

The angular velocity and period of rotation of a satellite (the period Is the 
time required to complete one revolution or 2 rr radians) are related by the for- 
mulas 

2iL. 

p “ p 


Where 


P = 


2n 


CJp 


P Is the period of the earth’s rotation, in dimensions of time; 
cop is the angular velocity of the earth to a processing equinox, 
obtained from the Appendix. 


Exercise A.7.2. Find the period of the earth’s rotation (a) in seconds, (b) 
in hours and decimal fractions of hours, (c) In hours, n)tnutes, and 
seconds. Compare with answers obtained in previous exercises. 


Astronomers call the average angular velocity of a satellite in orbit the mean 
motion. The formulas relating the period and th^ mean motion are 

r? _ 2w..„ 

P 



where 

r? is the mean motion In radians per second, or degrees per 
minute, or similar dimensions of angle per unit time; 

P Is the period in units of time. 


Exercise A.7.3. 

a) Find the mean motion of the orbiter when it is in a circular 
orbit of 160 n.mi altitude, in units of degrees per minute (.see 
exercise A.5.I.), 

b) Find the mean motion of the communication satellite of ex- 
ercise A.5.2, 
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The average angular position In the orbit Is called the mearr anomaly by 
astronomers. It Is found from the time and the mean motion by means of the 
formula 

M = f)(l'-t^)+ Ma 

Here 

M is the mean anomaly, in degrees, 

V is the mean motion. In degrees per unit time, 
t is the current time. In consistent units, 
to Is the reference time, called the epoch. 

Mo Is the mean anomaly at epoch. 

Exercise A.7.4. Suppose epoch is at midnight on July 18, 1981, and the 
mean anomaly at epoch is 40*. Find the mean anomally one hour and 
15 minutes after midnight for the orbiter when It Is In a circular orbit of 
160 n.ml altitute. 


AJ. MISCELLANEOUS PRACTICE EXERCISES 

Here are a set of practice exercises for which the formula to be used 
may not Immediately proceed the exercise. Vl/here not given, you are to 
find from previous pages the appropriate formula. Follow carefully' the 
last two calculator exercises. 


Exercise A.8.1. The orbital elements of the sun {pretend that the 
sun is In orbit about the earth) in the plane of its motion about the 
earth are: 


a -- 1.00000023 A.U, (Astronomical Units), 

c : 0.01671821 (this is eccentricity). 

(a) Find the semimajor axis of the sun’s orbit in n.ml. 

(b) Find the period of the sun (see u sun in Appendixi. 

(c) Find the mean motion of the sun In degrees per day. 


Exercise A.8.2. Ionized Hydrogen has a l,p of about 6000 secs. Find 
the value of the exhaust velocity In ft./sec. If an ionized hydrogen 
rocket motor developes 100 lbs of thrust, what is the rate of use of 
hydrogen, in Ibs/sec.? 
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Exercise A.8.3. Steam (H,0) has an l,p of 600 to 1000 secs, depend- 
ing upon the temperature. If steam at 750 *R develops an l,p of 800 
secs, how many lbs of water is used per siecond by a rocket motor 
developing 100 lbs thrust? What Is the effective exhaust velocity? 


Exercise A.8.4. The F y” function can be used to find the natural 
number e (this is not eccentricity) defined by 

e = J1 when x is very large 

Use this to find the limiting value of e, using x = 100, X = 1000, x 
= 10,000, X = 100,000. Compare the answers for each value of x. 


Exercise A.8.5. The square root of a number can be found by a cut - 
and • try method by means of the formula 

y _ Xn N 

" " ' 2X„ 

where 

N Is the number, 

X„ is the current trial value, 

X„ + 1 is the new trial value. 

Find the square root of 196. Suppose the first guess is 8. 

8 

STO 1 
2 

fy* 

196 

+ 

RCL 1 

-i- 

2 

-f- 


ANS 
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If this guess differs from the previous guess, use It to repeat the 
process. 

Exercise A.8.6. Write the calculator program for the cube root of 
81, if the formula Is 


y 2Xn + N 

^ “ 3x» 

and the symbols have the same meaning as given In the previous 
exercise. 


A.9. CAICULATOR FUNCTIONS 

For our purposes now, a function is a rule for converting one number (the 
number in the x-register, for example) to another number (which then appears 
in the x-register). Read Section 5 in the Owner’s Handbook, pages 85-107. A 
function may simply be considered a “more complicated” operator. The squar- 
ing operator is a function, since it is a rule for converting a number for another 
(for example, it is the rule for converting 3 to 9, 4 to 16, and so on). While we 
shall assume that the rule is programmed Into the calculator, and that it is ex- 
ercised by pushing the right buttons, we need to worry about two things; 

1. Does the number in the x-regIster have the right units (e.g., 
radians or degrees)? 

2. Does the rule “work” for the number in the x-regIster or 
y-register, in the case of two number functions)? 

Both concerns require that we think about the problem while performing the 
calculations. A good rule is to decide, before starting a calculation, on the 
dimensions to be used, Ihen convert all numbers during calculations. If there 
is a problem with this now, make all conversions before and write down the list 
of new numbers with their units. 

The calculator has 5 “named” functions. These are 

log, in, sin, cos, tan. 

The log and in functions require that the number in the x-register be a pure 
number (the number can be units of radians, but not degrees, and it can be in 
no combination of length, time, and/or mass.) This should cause no trouble, 
since we assume that the formulae have been stated correctly. The number 
appearing in the x-register cannot be negative or zero. As examples, use the 
calculator to find values for: 
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log 1, log 0, log - 1, log - 8 
1n 1, 1n 0, 1n - 1, in - 0. 

Note that whe 1 the rule does not “work", the calculator will announce your 
mistake. Other .obvious examples of “illegal" numbers are; division when the 
divisor Is zero, taking f VlTwhen the x-register is negative, and so on. They are 
“illegal" because the function is “not defined" in the mathematical sense for 
real numbers. 


The function “In" is called the “natural logarithm". It has as Its base the 
natural number e (Exercise A.8.4.). Use of calculators and computers makes the 
logarithm to the base 10, "log," far less useful than before. 


The sin, cos, and tan functions require that the number in the x-regIster be in 
units of degrees or radians. The calculator has a mode button h DEG | or h 


RADI which tells it what to expect (see page 93 of the Owner’s Handbook). 


To convert X degrees to radians, multiply the degree measure by 77/I8O. To con- 
vert radians to degrees, multiply the radian measure by 180 .. To convert an 
angle s pecified in degrees to radians on the calculator, key in the angle and 
press g |->R | . To convert a n angle specified in radians to decimal degrees, 
key in the angle and press f D-»|. 


Exercise A.9.1. 

a) Find the angle in degrees, minutes, and secs corresponding 
to 1 radian. 

b) Find the number of radians in 91 “ 40’ 23 ’. 

c) Convert the foliowing angles to degrees, minutes, and 
seconds: 


32.5042 degrees 
57.2175 degrees 
.01745329 radians 
1/9 radian 


Exercise A.9.2. Construct a tabie giving values of sin, cos, and tan for; 
o”, 15°, 30°, 60°, 75°, 90°, 105°, 120°, 150°, 180°, 

-0°, -15°, -30°, -60”, -75°, -90°, -105°, -120°, -150°, -180° 
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The radial distance of a satellite from the center of the earth Is given by the 
formula 

r = a <1 - e* ) 

1 + e cos e 

where 

r Is the distance from the center of the earth, In units of 
length, 

a Is the semimajor axis of the orbit. In units of length, 

e is the eccentricity of the orbit (dimensionless), 

© Is the true anomaly, the angle in the orbit between the 
satellite and the point of closest approach to the earth 
(perigee). 

Note that the altitude is given by 


H = r-R«, 


where 

H Is the altitude, 

is the equatorial radius. 

Exercise A.9.3. If an upper stage is in an elliptical orbit with a semima- 
jor axis of 13,185.39 n.mi and an eccentricity of 0.7267, what is the 
altitude of the satellite when e = 0“, 15°, 90°, 180°, 270°, 360°? 


Exercise A.9.4. The point of closest approach in an elliptical orbit is 
called "perigee," and the farthest point Is called the “apogee." Clearly, 
the denominator (1 -f e cos © ) must be the maximum, to make r 
smallest (at perigee), and smallest, to make r largest (apogee). At what 
value of © does perigee and apogee occur? Note that (1 e cos © ) 
is maximum at the maximum value of cos © . 

• 

The velocity of a satellite is given by the formula 

^ = xpTF]) 
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where 


r Is the radial distance from the center of the earth, units of 
length; 

a is the semimajor axis in units of length; 

n Is the gravitational constant, in units consistant with r, a, 
and the time. 


Exercise A.9.S. Find the velocity for the upper stage in the previous two 
exercises for values of © = 0°, 15°, 90°, 180°, 270°, 360°. What is the 
velocity at perigree? At apogee? Give dimensions of your answer. 


A.10. INVERSE FUNCTIONS 

A function Is a rule for converting one number to another. The Inverse func- 
tion is the rule for converting the second number back to the first. Another 
way of thinking: a function “does" and the Inverse “undoes." The function ap- 
pears in gold on your calculator. When it has an inverse, the inverse appears 
directly to the right of it in blue lettering. 

Most functions have inverses. A few do not. Examples of cases when the in- 
verse fails are: 

(a) division is the inverse of multiplication, except when the 
number is zero; 

(b) the f VT is the inverse of g x*, except when the number is 
negative. On the other hand g x* is always the inverse of f 
VX (Can you determine why?) 

The inverse of the In function is e\ As examples of its use do the following. 


3 

5 

7 

f In 

f In 

f In 

ge* 

ge“ 

ge* 


ANS 


ANS 


ANS 


3 


5 


1 


g e 


% 


96 


ge 


M 


f In 
ANS 


f In 
ANS I 


( In 


For a horizontal burn (that Is, gravity acts at right angle (90”) to or perpen- 
dicular to the direction of thrust), the weight of fuel burned In Increasing the 
speed of a rocket is given by the formula 

AV 

Wtb = Wp(e ^ - 1) 


where 

WpB Is the total weight of fuel burned, in lbs; 

Wp is the final weight of rocket and payload, after fuel has been 
burned, in lbs.; 

Av is the Increase* in velocity. In feet per second, av is the greek 
symbol pronounced "delta" and usually means "change in" 

or “difference " 

» 

c is the effective exhaust velocity of the burning fuel, feet per 
second. 


Exercise A.10.1. The o^biter has about 45 minutes to get far enough 
away from an upper stage so that the exhaust particles from the burn- 
ing rocket will not damage windows or tiles. To move 30 nautical miles 
away in 45 minutes requires that the orbiter increase its speed relative 
to the upper stage by 70 feet per second. If the Igp of the OMS fuel 
(Hydrazine) is 260 secs, and the orbiter weighs 200,000 lbs after the 
OMS burn, how much hydrazine is required, in lbs? 

ANS: 1680,8 lbs 

Exercise A.10.2. The OMS engines and the PCS thrusters both use 
Hydrazine, which has an l^p of 260 secs. The OMS engines develop 
12,000 lbs thrust total, and the PCS jets develop about 900 lbs. thrust 
per jet. Two jets can be used to separate from the Upper Stage. Which 
saves the most fuel, tt\e OMS or the PCS, for a Av of 70 ft/sec& (ANS: 
Both exactly the same. Explain.) 
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Exercise A.10.3 How many lbs. of OMS fuel is burned for each ft/ sec of 
Am acquired, when the orbiter weighs 200,000 lbs? When the orbiter 
weighs 160,000 lbs? 

• 

The amount of fuel required to perform a maneuver is important, since this 
fuel must be loaded before launch, and there ,is a limit to how much can be 
loaded, The flight planners must determine if enough fuel is aboard to perform 
the given maneuvers. 

The inverse problem is to find the Av when the weight of fuel burned is 
known. This is given by the formula 

AV = + c In (1 + 

Wp 

where the symbols have the same meaning as before. 


Exercise A.10.4. Find the Av when 1200 lbs of Hydrazine is burned dur- 
ing a separation maneuver, when the orbiter weighs 200,000 lbs. Note 
that In all of the above exercises, the burns were made in a horizontal 
direction. If the thrust acts against (or with) the pull of gravity, the for- 
mulas for AV do not apply. 




The AV required to change from one circular orbit to another close by is given 
by 

v.,„ = ViT a,* 

wnere ^ 

a, is the semimajor axis of the first circular orbit, 

a. Is the semimajor axis of the second circular orbit, 

M Is the gravitational constant, in consistent units. 

This formula applies when (Sj - a,) is less than 100 n.mi. 

Exercise A.10.5. The orbit'Hir is in a 160 n.mi. orbit and a payload is in a 
separate circular orbit. The orbiter is 10 n.mi. below the payload, that is, 

(a, - a,) = 10 n.mi. 
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a) How much av is required to transfer the orbiter into the 
same orbit as the payload? 

b) How many lbs of Hydrazine is required? 

c) Suppose the orbiter Is In the higher orbit, how does this 
change the answers? 


Exercise A.10.6. If the orbiter is in the same circular orbit as a payload, 
but several miles behind, what maneuvers would the orbiter make to 
catch up? 


A.11.THE TRIGONOMETRIC FUNCTIONS AND THEIR INVERSES 

The sin, cos, tan functions appear in gold letters and their inverses sin*', 
COS”’, tan”’ appear in blue letters to the right. As examples, run the following h 
DEG mode: 


15 

30 


60 

f sin 

f cos 


f tan 

g sin ”’ 

g cos 

-1 

g tan ”’ 

ANS 


ANS 


ANS 


The distance between two points on the surface of the earth is given by the 
formula (See page 97 of the Owner’s Handbook) 


where 


D = 60 cos ’ (cos (lat,) cos (lat,) cos (long, - long,) 
+ sin (lat,) sin (lat,) ) 


D is the distance, in nautical miles; 
lat,, lat, are the latitudes of the first and second points respec- 
tively, with southern latitudes entered as negative 
numbers; 

long,, long, are the longitudes of the first and second points respec- 
tively, with eastern longitudes entered as negative 
numbers. 


29 





ExerelM A.11.1. Two tracking stations, one near Los Angeies and the 
other near Baytown, Texas have the locations; 


Los Angeles 

34” 20’ 42” N 

lat 


118” 42’ 05” W 

long 

Baytown 

30*41’ 05” N 

lat 


95* 06’ 32” W 

long 


What is the distance between the tracking stations? 

In solving for the location of an object in orbit, Kepler Introduced an Angle E, 
called the eccentric anomaly. The formulas relating the eccentric anomaly to 
the true anomaly are 


E = 2 tan-' 



0=2 tan- 


-1^-^ tan EJ2 1 
1 - e 


Exercise A.11.2. For the orbit in exercise A. 10.5., find values of E 
for each value of E. How much do E and 0 differ at perigee? At 
apogee? At 0 = 90°? 


A.12. QUADRANT CHECKS 

If the circle is divided into four equiangular parts 

Part I O'-QO" 

Part II gOMSO” 

Part III 180“-270" 

Part IV 270 ”-360“ 

each part is called a quadrant. Quadrants are r^ombered counterclockwise 
from I to IV. 

Exercise A.12.1. With the above system of numbering quadrants, what 
is the logically consistent way to label the direction of increasing 
angles? Decreasing angles? 
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One problem with trigonometric functions Is the Inability to find the angle us- 
ing the Inverse function when the angle lies in the second or third quadrant. 


Examples are h Deg 
46“ 

node: 

135“ 

135 “ 

225* 

sin 

sin 

cos 

cos 

sin-’ 

sin-’ 

cos’ 

cos-’ 

1 ANS 

1 1 

IansI 

Ians! 


The situation is that the inverse function aiways yieids an angle called the prln- 
cipal value. 

I nverse Function Principal Value 

sin-’ -90“ to -j-90" 

cos ’ 0 to 180 

tan-’ -90 to -i- 90“ 

If we have only the value of the sum of the angle, or if we have only the value of 
the cosine of the angle, we are stuck with the principal value. 

Exercise A.12.2. Using the tables constructed in Exercise A.9,2., find 
corresponding inverses and compare with the original angles. 

If we have both the sine value and the cosine value of the angle, we can find 
the value of the angle In the Interval 0 to 360“ by using the inverse tangent 
function and the truth table. 

Given: sin e , cos o 

Find: e = tan-’ 

cos 6 


COS>S^ 

-t- 

- 

+ 

0“ 

-♦■360“ 

- 

-M80“ 

-♦-180“ 
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The truth table is used as follows, 


1 . 

2 . 

3. 

4. 


If sin e Is positive, and cos 0 
If sin e • is positive, and cos e 
If sin 0 Is nbyatlve, and cos 0 
If sin Q is negative, and cos o 


is positive do nothing to e , 
is negative, add 180* to 9 . 
Is positive, add 360* to 0 . 
is negative, add 1^* to 9 . 


Exercise A.12.3. Using tables constructed In Exercise A.9.2., use the 
truth table to find the angles from the Inverse tangent. 


Another problem with trigonometric Inverses Is that we can never discover If 
the original angle was outside the internal 0 to 360*. As examples, show that 
the truth tables yield the same result for the following angles 

9 = 135, 495°, 855*, 

9 = -225,-585". 

That Is, taken sin © , cos 0 , and find 0 = tan-' sin 0 and apply truth table. 

cos 0 

There Is apparently nothing we can do about this, so we adapt to it. 


A.13. EVALUATION 
A.13.1. Learning Objectives 

The principal objectives of Section A were threefold: first, to teach basic 
computional skills through use of the calculator; second, to explore calculator 
capability and limitations in examining formulas and functions used in flight 
planning; and thirdly, to teach the basics of problem analysis, relating word 
sets to formulas to calculator logic. 

Having completed this section and related references, you should have an 
understanding of: 

• the logic and capability of the calculator In solving simple 
problems Involving basic mathematical functions 

• the basic components of formulas and the rules for problem 
solving 
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• the calculator functions which simplify problem solving 

• basic mathematical and space flight related terminology. 


A.13.2. Outcome Measures 

You should now be able to; 


1. Write the calculator code or key steps In evaluation of a for- 
mula; 

2. Read a simple problemi statement and select given variables 
and those to be calculated; 

3. Evaluate formulas with the calculator containing powers, 
roots, logarithms, exponentials, trignometric functions, or In- 
verse functions, or all; 

4. Make use of calculator conversions, such as radians to de- 
grees, fractions of degrees to minutes, etc,; 

5. Make conversions of units of measure from one set of mea- 
surements through use of equivalent relations; 

6. Demonstrate a passing acquaintance (they will become 
clearer as we proceed) with the following terms. 


radius 

velocity 

angle 

degree 

revolution 

rocket 

conversion 

formula 

mean motion 

mean anomaly 

elliptical orbit 

astronomical units 

nautical miles 

specific impulse 


altitude 

force 

radian 

period 

satellite 

thrust 

dimension 

equation 

transfer 

true anomaly 

eccentricity 

semimajor axis 

angular velocity 

gravitational constant 


axis 

mass 

rotation 

epoch 

orbit 

launch 

function 

impulse 

payload 

perigee 

eccentric anomaly 
international feet 
exhaust velocity 
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B. PROGRAMMING THE CALCUilATOR 


BJO. INTROOUCTSON 

You have been solving problents manually on the calculator using a step-by- 
step process. In this section you will begin to write programs for the calculator 
that will enable you to press a single key on the calculator, commanding it to 
execute all the steps necessary for the solution of a problem. It is not quite 
that simple, however. 

But for you now a program will consist simply of a series of calculator 
Keystrokes pressed to solve a problem manually. The calculator will rerriember 
these keystrokes when keyed in, then execute them in given order at the press 
of a single key. We introduced use of the storage key in Exercise A.8.5. Pro- 
grams are extremely useful when solving problems involving use of the same 
process (or formulas) but containing different vaiues for the variables. 

You must learn how to construct programs, load them, run them on the 
calculator, and record them for future use. Some programs are simple to write, 
others are not. Understanding your calculator will also serve to provide the 
foundation for future work in programming larger computers. 


B.1. PROQRAMMINQ PROCEDURE THE JULIAN DAY NUMBER 

Review the material on pages 17-23 In the HP-67 Owner’s Handbook and Pro- 
gramming Guide, and do the suggested exercise. This material provides an in- 
troduction to running a prerecorded program, creating, loading, and running a 
program, and recording of the program. Now turn to Section 6, pages 123-145, 
entitled Simple Programming, and do the suggested activities. Spend some 
time learning the calculator code, pages 129-130, 

A stored program is relatively simple. It consists of code 

g LBL f (I) or f LBL (L), 
set of instructions, 
h RTN. 

The first entry is a LABEL designation, and the last entry is a RETURN (RTN). 
The first identifies the beginning of program or first step of program memory 
through the label L (one of letter keys A through E) or label I (one c ' letter keys 
a through e), and the first RTN identifies the end of that program. These are re- 
quired by the calculator control logic, otherwise the calculator could not begin 
or terminate. 


34 



To illustrate programming construction, we will write a program for com- 
puting Juilan Day Numbers from the Calendar Date. Note that writing a stored 
program Is a form of problem-solving since we must prepare: 

(1) a list of data that will be contained In the answer, together 
with dimensions; 

(2) a list of input, or "givens," together with dimensions, from 
which the answer must be constructed; 

(3) a list of formulas connecting the given data and the 
answers, or output. 

The plan of a solution, in programming construction, is called a "flow 
diagram" or flow chart. Such diagrams are useful In solving programming 
problems, but their use is not restricted to this area. It is preferable to some to 
work out the calculation process, constructing a flow diagram as they go, and 
then write the program. Others draw a flow diagram after the program Is com- 
pleted for reference purposes. The most important standard symbols used are: 


for input and output 
or 

general operations 



for decisions 


These symbols are joined by lines, with logical flow normally from top to bot- 
tom unless indicated otherwise by arrows. Flow diagrams are discussed on 
pages 140-145 of the Owner’s Har?dbook. 

A "subroutine" is a section of a program that is performed as a separate en- 
tity, and stored ready for use when needed in a master program. It Is often con- 
venient to write a very complicated program as a set of subroutines or sub- 
programs. For our purposes a subroutine will be no more than evaluation of a 
single formula. 
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While a stored program Is simple, the preparation tor and documentation of 
a program to be stored may not be simple. Unless certain rules are strictly fol- 
lowed, chaos may result. The procedure Is as follows; 

1. Define; • Input data 

• Output data 

2. Select dimensions for • Input data 

• Calculations 

• Output data 

3. Write down formulas. 

4. Note logic tests required, branches, and subroutines re- 
quired, If any. 

5. Do flow diagram of calculations. 

6. Test the code through; • Editing 

• Using known results from other 
progran«s. 

• Exerckuing all paths 

7. Document the program by; • Title 

• Date written 

• Brief description of what 
program does 

• Formulas used 

• Operating limits and warnings 

• Flow diagram 

• Operating instructions 

• Code and notations 

For astronomical considerations time is an important factor. We can 
describe the positions of celestial objects in terms of time; we can describe 
<he positions of earth orbiting bodies in terms of time. A system of Julian Day 
Numbers was introduced to reduce computational labor in the comparison of 
time spans made many years apart and to avoid time ambiguities. The Julian 
Date, which we will use in other calculations, is given for every day of the year 
in the American Ephemeris and Nautical Alamanac. Its calculation, however, 
will be helpful. The Julian calendar epoch was chosen to be January 1, 4713 
BC, and it assumes the length of a year to be exactly 365.25 days. The Julian 
Date may be calculated from the relation 

JD = (365.25) y + (30.60) m -f Day -t- 1720981.5, 
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where 

y = year -1, if month is less than or equal to 2 ( < 2), 

y = year, If month is greater than 2 ( > 2), 

m = month + 13, If month Is Jess than or equal to 2 ( s 2), 

m = month + l. If month Is greater than 2 ( >2 ). 

This relation has conditional statements; depending on the month, the calcula- 
tion may go one of two ways, with variables y and m given by 

y = year -1 
m = month 13, 
or 

y = year, 
m = month + 1. 


In a flow diagram, this would be described by 

day 

month 

year 


decision 



m = month 13 


y,m 

JD = JD (day, month, year) 
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Let us now construct the program documentation. With a present caiendar 
date (day, month, year) as input and Julian Date as output, we start with the 
preliminaries, program description and formula. Follow carefully this 
documentation. Key in the program as coded. 
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Prosnim Doeumantatlon (1) 

(Title, Date, Brief Description, Formulas, Limits) 

Title: Julian Day Number from Calendar Date 

Date: 

Name: 


Program Description: 

Program computes Julian Day Number for 00 Mrs. 00 Min. 00 Sec foi In- 
put calendar date 


Equations, Variables, etc: 

Julian Day No. = (365.25) y -»■ ,(30.6001)m + Day 1720981.5 

where 

y = year - 1, If month less than or equal to 2 ( s 2), 
y = year, If month greater than 2 ( > 2) 
m = month, + 13, If month less than or equal to 2 ( s 2) 

M = month + 1, If month greater than 2 ( > 2) 


Operating Limits: 

Valid for calendar dates March 1, 1900 to February 28, 2100 
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Program DocumenfafSon (9 


(Flow Diagram) 



RTN 

I 

I 
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Prosram Doeumwitation (3) 
(Operating Instructions) 
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Program Documentation (4) 
(Code and Notations) 
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Exercise B.1.1. Write and document below a subroutine for calculating 
the altitude of a satellite from its orbital parameters and the true 
anomaly (Exercise A.9.3). Use documentation sheets which follow. 


Exercise B.1.2. Write and document subroutines for calculating: 

a) the velocity of a satellite from its Oibital parameters and the 
true anomaly (Exercise A.9.5.), 

b) the velocity at perigee and the velocity at apogee from orbital 
parameters (Exercise A.9.4.), 

c) the period of a satellite (Exercise A.5.1.), 

d) the mean motion (Exercise A.8.1.), 

e) the mean anomaly from the mean motion (Exercise A.7.4.). 
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Program Documantatlon (1) 

Title: Altitude of Satellite from Orbital Parameters 

Date: 

Name: 


Program Description: 


Equations, Varlabies, etc: 


Operating Limits: 
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Program Domiintntatton 

TItta: Altitude of Satellite from Orbital Parameters 


I 

a 


I 


I 

a, e, 0 , R«, 


INPUT 




Prognm DocumwtoMen @ 

THI*: Altitude ot Satellite from Orbital Parameters 



46 














47 










POS^N IN PLANE OF ORBIT FOR NEARLY CIRCUUR ORBITS 

For small eccentricities (e), the true anomaly maybe found from the formula, 

e = M 2e sin M -I- 0 * 8ln 2M 
+ (13 sin 3M - 3 sin M) 

The position vector of a satellite can be found at any time t by the following 
procedure. 

Given; Mg, tg, a, p , e and time t, 

Compute: 

O tr a®'* 

p = (Exercise A.5.1.) 

"vT 

r? ss (Exercise A.8,2.) 

M = ^ (t - tg) + Mg (Exercise A.8.3.) 

Test: 

If e > 0.1 STOP with e In X-Reg. 
e < 0.1 GTO next step. 


Compute: 0 from e and M using formula above. 
Compute; 


r (Exercise A.9.3.) 

1 + e cos 0 ' ' 


V = 



— (Exercise A.9.5.) 


The output data: P, 6 , r, V at time t. 


Exercise B.2.1. Write and document a calculator stored program for ob- 
taining 0 , r, and V at time t. (Specify dimensions of output!) 
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a and M are input from Input data 
for calculation of P 


calculated P is input for calculation of n 


calculated v Is input for calculation of M 
-t, to and Mo from Input data 


calculated M is Input for calculation of 
0 - e from input data 


calculated 0 is input for calculation of r 


calculated r is input for calculation of v 













Exercise B.Z2. Modify the above program so that the output is the 
altitude H, as weii as 0 , r and V at time t. 


Measures of a distance r and an angie 0 are caiied the "poiar coordinates” 
of a point. Topics invoiving motion about a point can better be investigated 
with poiar coordinates rather than by rectangular coordinates (ordered pairs of 
real numbers (x, y)). Graphs of polar coordinates are generally done on polar 
graph paper, paper with angular and radial divisions. To make a graph we mark 
the point for each order pair (r, 0 ). For example, for r = 1 + cos 0 , a table 
could be constructed from rwalues calculated by inputting different values of 
©. 


0 

0° 

45° 

90° 

120° 

150° 

180° 

210° 

240° 

270° 

300° 

330° 

r 

2.0 

1.7 

1.0 

0.5 

0.1 

0 

0.1 

a5 

1.0 

1.5 

1.9 


Using polar graph paper, find the point for each pair (r, 0 ), see if a pattern 
develops, and then join the plotted points with a smooth curve. 

Exercise B.2.3. Using the calculator program from the previous exer- 
cise, lei 


to = 0 
Mo = 0 

Find r and © for t = 0. Input t in 15 minute intervals, draw up a table, 
and plot the values of r vs 0 on polar graph paper for: 

(a) the orbiter in a 160 n.mi. circular orbit (e = 0.0), 

(b) the orbiter in an orbit with a = 3875.1944 n. mi. and e = 

0.07, 

(c) the orbiter in an orbit with a = 4004.3676 n.mi. and e = 0.1. 

Plot all on the same graph paper. The graph will show a circular orbit 
and two elliptical orbits with different eccentricities. 

Exercise B.2.4. Find the velocity for each orbit above when all are coin- 
cident (take program V output at 0 point seen in graph). The dif- 
ference in velocity is the av required to transfer from one orbit to the 
other at the point of coincidence. 
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Exercise B.2.5. Write a program to calculate amount of fuel in pounds 
when the 1^,, A v, and orbiter weight are given. How many pounds of 
hydrazine is required to transfer between any two orbits above if the or- 
biter weighs 200,000 lbs? (See Exercises A.6.2. and A.10.1.) 


8.3. ELMPTICAL ORBITS 

In Exercise B.2.3 a set of orbits was plotted. The first orbit was circular, but 
the more general path followed by an orbiting body is elliptical. By 
mathematical definition an ellipse represents the set of all points (x, y) in an x, 
y plane, the sum of whose distances from two distinct fixed points (foci) is con- 
stant. If we fasten the end of a fixed line of string to thumb tacks and draw the 
string taut with a pencil, then the traceable path of the pencil will be an ellipse. 
The thumb tacks in this case would be the focal points (foci) of the ellipse. 

In Figure 1B(a), C is the center of the ellipse and points A and B are the ver- 
tices of the ellipse. The points P and Q are called the focus points, or the foci 
of the ellipse. In every ellipse, the foci are located along the major axis be- 
tween the vertices and the center. 

The length of the semimajor axis is generally denoted by the letter a as 
shown in Figure 18(b). Thus the length of the entire major axis is 2a. The 
length of the semi -minor axis is generally denoted by the letter b. The length 
of the segment from the center to one foci is generally denoted by the letter c; 
thus the length of the segment from one foci through the center to the other 
foci is 2c. 

Since the foci are located along the major axis between the vertices and the 
center, c is grfe-jter than zero but less than a (0 < c < a), and therefore as c 
becomes small approaching zero, the foci approach the center and the ellipse 
becomes more circular. As c approaches a in length, the foci approach the ver- 
tices and the ellipse flattens out. Thus, the shape of an ellipse can be de- 
scribed in terms of the relative sizes of a and c. This common measure of the 
shape of an ellipse is called its eccentricity, given by the ratio of c to a. 

Six orbital elements are used in describing a satellite’s path in space. Three 
of them define the orientation of the orbit with respect to a set of axes, two of 
them define the size and shape of the orbit, and the sixth, with the time, 
defines the position of the body within the orbit at that time. For a satellite of 
the earth traveling in an elliptic orbit, the fundamental reference plane is the 
equator. 
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The six elements defining the orbft and the position of the body in It are the 
semi • major axis and eccentricity of the elliptic orbit, longitude of the ascend- 
ing node (to be discussed in later sections), argument (angle) of perigee, in- 
clination of the orbital plane toward the equator, and the time of the passage 
through perigee, or the mean anomaly at the initial moment of time. Perigee, 
apogee, and true anomaly are described In Figure 2B; argument of perigee and 
angle of inclination are described in Figure 3B. 



Figure 2B. An elliptical orbit. For a satellite at position S, r is the radius vector 
from the center of the earth Q. P, the point of closest approach to earth, is 
perigee. A is apogee. PA is the iine of apsides. The angle © , measured be- 
tween QP and QS, is called the true anomaiy. 
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Figure 3B. Earth orbit. With P denoting perigee, g is the argument of perigee 
for an orbit with angle of Inclination I. 


B.4. KEPLER’S EQUATION 

Johannes Kepler (1571-1630) formulated three major laws of planetary mo- 
tion. The first was that the orbit of each planet is an ellipse with the sun 
always located at one focus. But a fundamental problem of motion for a 
satellite in any elliptical orbit (earth about sun or object about earth) Is to find 
its position at any given time. Kepler worked on this problem for twenty years, 
then gave up. 
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The solution that he found and thought unsatisfactory, however, has not 
been improved upon In the more than 300 years since his death. It was derived 
from his second law which Indicates how the speed of a satellite changes in 
its orbit. Kepler defined a new angle, the eccentric anomaly, by means of an 
auxilliary circle of radius equal to the semimajor axis and a geometric con- 
struction, as shown In Figure 4B. 



Figure 4B. Eccentric Anomaly. Constructing a circle about the elliptical center 
C as shown, the eccentric anomaly E can be defined through the true anomaly 
0 . 
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The formula derived is 


M = E - e sin E 


vt/here 

M is the mean anomaly, 

E is the eccentric anomaly, 
e is the eccentricity. 

If M is known, E may then be found by approximation methods. We wiii use 
that which is called the Newton-Raphson method. We use the relation 

p _ p j. M - E„ -t- e sin En 
+ 1 - e cos E„ 

where n is a number starting with zero representing the n"’ approximation. The 
o'^ approximation for E is M for eccentricities smaller than 0.1. That is to say, 
fore < 0.1, let E« = M. 

The first approximation for E is 

_ M-Eo + e sin E^ 

’ “ ^ ^ 1 - e cos Eq 

If the difference E, - E„ is large, we take the second approximation for E: 

___ M-E, + e sin E, 

* ~ ^ 1 " e cos E, 

We then check the difference Ej - E,. If that is large we continue with E3, etc., 
until the difference E„ + , - E„ is so small that further calculation is un- 
necessary. A calculation such as this, while appearing complex and lengthy, 
may easily be done by the calculator through programming. We will do this in 
the next exercise. 

Having obtained a value for E, the true anomaly can be calculated from 
with 

V = 
and 

M-= n (t - to) -H M„ 
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Kepler’s equation is solved from a given t, t,, Mo and v . Programming for a 
solution Is given in the following exercise. Since it is a more complex exercise, 
the steps will be provided. Study the flow diagram, and give attention to the 
approximation or iteration procedure for determining E. 


Exercise B.4.1. Write a program to solve Kepler's equation. 


Inputs: t, to minutes 

Mo degrees 
a n.mi 

Constants: ^ (n.mi)*/sec’ 

60 sec/min 

180 for deg to radian 


Computations: 


T? = (60) rad/min 

M = n (f - to) + Mo ( 7T /180) rad 


E 


n + 1 


M - Eq e Sin En 
^ 1 - cos E„ 


0=2 tan-’ tan E/2) deg 


Storage: M, e, E„, e 

Working storage required: 4 or more 


Tests: test on difference E„ + , - E„ and stop if difference less than 10-* 
radians 
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rtoonim Doeumentatitm (1) 

(Title, Date, Brief Description, Formulas, Limits) 


TKIe: Kepler’s Equation 

Date: 

Name: 


Program Description: Program solves Keplers equation for Mg eccentric 

anomaly to obtain the true anomaly from given input data, t, L, Mg, a, and M 


Equations, Variables, etc: 


M = E » e sin E 


=^/i 


M = TJ (t - g 4- Mg 


" 1 - e cos E„ 


0 = 2 tan"’ / fjjt-- 

aV'-®" 2 


Operating Limits and Warnings: 


Use M for Eg in En + 1 = E„ 4- 


M - En 4- e sin Eg 
1 - e cos E„ 


Stop calculation for E„ ^ , when E„ + , - E„ < 10" 


M = T? t 
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Program Documantation 

(Title, Date, Brief Description, Formulas, Limits) 

Title: Kepler's Equation 
















Program Pocumontathm (£) 


Tlllo: Kepler’s Equation 

































































Program Ooeumwitatton ( 4 ) 


Title: Kepler's Equation 


Title: Julian Day Number from Calendar Date 























































































ExwelM B.4.2. Compare the values of 9 obtained from the program of 
Exercise B.4.1. with values of 9 obtained from the program of Exercise 
B.2,3, for values of e = 0.05, e = 0.1, e s 0.2, e 0,5, 


Exercise B.4.3. Using the calculator program of Exercise B,4.1, plot 
values of r vs 9 on polar graph paper for the following oibits. 

(a) A circular orbit of 160 n.ml altitude (e = with i, = 0, M, s 
0 for t at intervals of IS min, 

(b) An elliptical orbit, with a = 131^,3858 n.ml, e = 0.7266723 
and to = 0, Mo - 0 at inten/als of 30 n.mi 

(c) A circular orbit, at an altitude of 19322.91 n.ml, and t„ = 0, = 

180* at intervals of 1 hour 


B.S. EVALUATION 


B.S.I. Learning Objectives 

The principal objective of Section B was to teach more complex problem solv- 
ing based on task analysis which involved examples of flight planning problems 
with arbitrary and incremental variable input. The complexity of the task was 
simplified through the use of the calculator programming capability. 

Having completed study of this section and related references, you should 
have an understanding of: 

• the conceptual content of problems 

• the types of data needed for problem evaluation 

• the use of multiple formula in problem evaluation 

• the logic of flowcharting 

• the methods of constructing a calculator program 
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B.5.2. Outcom* Meaiuic* 

You should now bo able to; 

1. Identify Input data and dimensions needed to solve prob- 
lems, 

2. Identify appropriate formulas for problem solution, 

3. Flowchart a problem to solution, 

4. Write a calculator program for problem solution, 

5. Document and edit programs, 

6. Demonstrate a familiarity with orbital terminology. 

If not, QTO B.1. 
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a POSITION AND VELOCITY 
IN THE ORBITAL PLANE 


C.O. INTRODUCTION 

The purpose of the following is to develop ways to calculate the position of a 
satellite in the plane of its orbit about the earth when the orbits are highly ec- 
centric as well as circular. The given forms of input data will be typical of that 
encountered for satellites or otnei orbiting vehicles in current mission plan- 
ning. Because of the importance of upper stage delivery in orbiter flight pro- 
grams, we will look at the details of orbit transfers in the plane. 

C.1. LIST OF SYMBOLS AND DIMENSIONS 

A listing now will be made of symbols used to denote all the variables to be 
employed in our calculations. Many of them have been introduced before. 
Here, dir tensions will be given using the three basic units of mass (M), length 
(L), and time (T). 

V = velocity, dimensions L/T 

r = distance from center of earth, L 
a = semimajor axis, L 
^J = gravitational constant, UIV 
e = eccentricity, dimensionless 
rp = perigee distance from center of earth, L 

r^ = apogee distance from center of earth, L 

H = altitude above equatorial radius, L 
r; altitude at apogee, L 
Hp = altitude at perigee, L 
P = period, T 

V = velocity in orbit, L/T 

Vdro = velocity in circular orbit, L/T 
Rpq = equatorial radius of earth, L 
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Vf, = radial component of orbital velocity, LIT 
= horizontal component of velocity, L/T 
p = parameter of the ellipse 

V = mean motion, L/T 

O - true anomaly 

E = eccentric anomaly 

M = mean anomaly 

7 = flight path angle relative to local horizontal 
1^, = specific impulse, T 

F= rocket thrust, M L/T 
Wb = rate of fuel burned, M LIV 

WjB = total fuel burned, M L/T 

Av = velocity change, L/T 

C = effective exhaust velocity, L/T 

C^. BASIC FORMULAS FOR THE ELLIPTIC ORBIT 

Again, while a number of the basic formula have been introduced, a listing of 
all formulae necessary for a description of an object in o'-bital flight is 
presented below. 

(1) g = (gravitational) 

^eq 

The legal value of g is 32,1740486. (Note: the legal value of g is used to convert 
weight in lbs to mass in slugs.) 

(2) P = 

(3) r? = 

( 4 ) r, ^ 


2 rr a' 


3/2 




2 n 
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(5) 

P = 

2 fr 
T? 

(6) 

H = 

r- R«, 

(7) 

r = 

H + 

(8) 

M = 

(t - g + M, 

Computation of Semimajor axis; 



(9) 

a = 

P 

1 - e’ 

(10) 

a = 

fa 

1 + e 

(11) 

a = 

1 - e 

(12) 

a = 

2r,-p 

(13) 

a = 

[e 

2fp - p 

(14) 

a = 

fa + fp 
2 

(15) 

a = 

r« + fp 
2 

(16) 

a = 

r 

2 - rV» 

(17) 

a = 

(PVI? ) 2/3 
2 
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Computation of eccentricity: 


(18) 

e = 

'a + fp 

(19) 

e = 


(20) 

e = 


(21) 

e = 

A-1 

a 

(22) 

e = 

a 

(23) 

e = 

1- P 

ra 

(24) 

e = 

p -1 

Computation of parameter p: 



(25) 

P = 

1 

CD 

(26) 

P = 


(27) 

P = 


(28) 

P = 

r£p 

a 

(29) 

P = 

re (1 - e) 
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(30) 

P = 

•"p (1 + e) 

(31) 

P = 

2r.rp 

r, + Tp 

(32) 

P = 

r*V^ r*v*cos’ y 

Computation of radial distance r; 



(33) 

r = 

a (1 - e») 

1 + e cos © 

(34) 

r = 

P 

1 + e cos 0 

(35) 

r 

2a 

I “ 

1 + aV» 

(36) 

- 

a(l± / 

V \ cos* 7 

(37) 

r = 

a (1 - e* ) tan y 

e sin 0 

(38) 

r = 

a (1 - e cos E) 

(39) 

r = 

V7T a e sin E 

(40) 

r = 

VTp Vh 


Computation of radial distances at apogee and perigee; 
(41) r, = a (1 + e) 
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(42) 

fp = 

a (1 - e) 

(43) 

r. = 

ap 

^P 

(44) 

r, == 

2a - fp 

(45) 

r. = 

P 

1 - e 

(46) 

r. = 


(47) 

r. = 

Q. 
CL 1 

CM 

(48) 

= 

?P 

(49) 

fp = 

2a - Tg 

(50) 

r — 

P 

'P 

1 + e 

(51) 

fp = 


Computation of true anomaly: 


f 

(52) 

0 = 

2 tan-^ ( / 


tan -|) 


(53) 

(54) 

(55) 


0 



0 = 


© = 




1 - sin E 
1 - e cos E 
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(56) 


= C 08 ’ (y (COS E - e) ) 


(57) 


0 = cos ' ( ^ ) 
' er ' 


Computation of eccentric anomaiy: 


(58) 

(59) 

(60) 

(61) 


E= 2tan-M ^/T^tanf) 


E = COS'* (1^1^) 

0 


c oin -1 / r sin 6 V 


E = cos- 


1 


aV^- 1 
JL_ 


® \a^ + 1 


(62) 



(1 -e») 
cos* y 


Computation of fiight path angie: 


(63) 

y = 

Sin-' (^) 

(64) 

y = 

cos ' — 

V 

(65) 

7 = 

Vh 

(66) 

7 = 


(67) 

7 = 

cos-M >f) 
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(68) 

7 = 

tan-’|’(1 - p) tan ©j 

(69) 

7 = 


Computation of velocity 



(70) 

V = 


(71) 

Vcl. = 


(72) 

V = 

I + B COS E) 

a (1 - e COS E) 

(73) 

V = 

I UP 
V r cos 7 

(74) 

V = 

> 

/VS + VS 

Computation of radial and horizontal velocities; 

(75) 

Vh = 

V cos y 

(76) 

V. = 

V sin 7 

(77) 

= 

(e sin 0 ) 

(78) 

Vh = 

>/r 


Escape velocity computation: 


(80) 


( 81 ) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 


I.P = F/Wp 

F = Wp l.p 

Wp = F/!^ 

I.p = c/g 

Wtb = Wp(e-^-l) 

Av = cln(1+^) 

^v,,rc =V7T 


Kepler's equation; 

(87) 

( 88 ) 


M = E - e sin E 

M - E„ + e sin E,, 


^ + 1 ~ Ffi -f 


1 - e cos En 


Series expansions for low eccentricities; 
(89) 




+ ( 
+ (■ 


5e» 

lle^ . 

17e® 

4 

24 

192 

13e» 

43e“ 

95e' 

12 

192 

512 

103e‘ 

451 e® 

4 . 

96 

480 

... 


- ... ) sin 2 M 

- ... ) sin 3 M 
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, /I097e* 5^7®' , 


+ - ...) sin 6 M 




(90) 


E = M + (e 


§1 . e° e^ 

8 192 9216 


+ ... ) sin M 


+ (f - | + |-...)sin2M 


1 / 36 * 276 * , 2436 ^ \ o Kji 

W + 5120 - ■>*3M 


+ + ... ) sin 4 M 

3 15 


+ ( _ 3125®^ V , 5,^ 

^ ' 384 9216 ^ ■ • •' ^ 


+ ■ .) sin 6M 


C.3. TRANSFER OF ORBIT 

All space vehicle trajectories, at least during some phase of the vehicle’s 
mission, can be classified as transfer orbits. This general classification is 
basic to most missions, because, In almost all instances after the space ve- 
hicle has been established in an initial orbit, a velocity impulse is given to the 
vehicle \which causes it to transfer to another desired orbit. Transfer orbits are 
employed because a direct launch into the final orbit is almost always 
physically Impossible and because an Intermediate orbit can be more 
desirable for initiating a precise transfer to a translunar or interplanetary orbit 
than a direct transfer. 


73 



In consideration of an orbital transfer It Is necessary that the orbital 
characteristics of both initial and final orbits be known, With that information 
it can be determined whether a transfer orbit is necessary and the method by 
which transfer is to be accomplished, 


Exercise C.3.1. Radar observations indicate that an orbiting satellite 
has an aiititude at apogee of 600 n.mi and a perigee altitude of 150 n.mi 

a) Find the semimajor axis (Eq, (15) ) 

b) Find the eccentricity Eq. (19) or (21) ) 

c) Find the period (Eq. (2) ) 


Exercise C.3.2. The orbiter is in an orbit with a semimajor axis of 
4454.91 n.mi and eccentricity of 0.2, 

a) What is the perigee altitude? 

(Eqs. (6) and (42)) 

b) What is the apogee altitude? 

(Eqs. (6) and (41) ) 

c) What is the radial velocity at perigee? at apogee? 

(Eqs. (70) and (76) ) 

d) What Is the horizontal velocity at perigee? at apogee? 

(Eqs, (70) and (75) ) 

e) What is the radial distance when the eccentric anomaly is 

90”? (Eq. (8)) 

f) What true anomaly coincides with this eccentric anomaly? 

Transfers orbits may be classified In several ways. One of the methods of 
classifying transfer orbits Is by the type of thrust used to either initiate or ac- 
complish the change of trajecton/. The transfer could be the result of a single 
,/elocity impulse, a sequence of velocity impulses, a low impulse thrust, a high 
Impulse thrust, or a continuous thrust along the transfer trajectory. Engines 
may be fired toward the earth, adding to the velocity in a direction away from 
the earth, they may be fired in the opposite direction, or in any direction in be- 
tween. 

Velocity is really a vector quantity. Some quantities, such as mass and time. 
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can be represented by a single real number scaled to an appropriate unit of 
measure. Such quantities are referred to as scalar quantities, and the real 
number associated with each Is a scalar. Some quantities can not be de- 
scribed fully by a single real number because they have both magnitude and 
direction. Force (thrust), and displacement are examples In addition to velo- 
city. The mathematical object used to describe each vector quantity Is called a 
vector. 

A vector Is represented geometrically by an arrow. By definition, a vector Is 
a directed line segment, with a tail end (initial point) and a head end (terminal 
point), and it Is understood to be directed from tail to head. The length of a 
vector represents Its magnitude and the direction of the arrow represents the 
direction in which the vector acts. 

On the earth, or from the earth surface reference, we make use of the up- 
down direction, which we designate as vertical. A direction at right angles to 
the vertical, perpendicular to the vertical or parallel to the earth surface. Is call- 
ed horizontal. Often, a vector quantity is described in terms of vertical and 
horizontal components as shown in Figure 10, 

If the earth surface is used as reference, a satellite In circular orbit would 
have a velocity that was always horizontal, perpendicular to the vertical which 
will be referred to as the radial direction. In an elliptical orbit, however, the 
velocity is completely horizontal only at apogee and perlgea. Throughout all 
other points of the elliptical orbit the satellite has radial and horizontal com- 
ponents as shown in Figure 2C. 


Exercise C.3.3. The orbiter is launched on a trajectory such that at the 
completion of the OMS 2 burn it Is traveling horizontally with velocity of 
25,495.83 feet/sec at an altitudn of 607,611.55 ft. 


a) 

What 

Is the semimajor axis? 

b) 

What 

is the eccentricity? 

c) 

What 

is the perigee distance? 

Hint: 

The 

data given is Vh, Vr and H. 


The 

data wanted is a, e, rp 


Use 

(74) to find V (Vr = 0 from above). 


Use 

(7) to find r. 


Use 

(16) to find a. 
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Vertical 

(radical) 

direction 


V 






Or 


(a) 


vector V 



Horizontal 

direction 






















Vh 


Horizontal 


(b) 

Figure 1C. Vector components, (a) The vector labeled V Is described by Its 
magnitude, V, and direction, a , from the horizontal, (b) The vector V is 
equivalently replaced by two vectors, a horizontal vector with magnitude 
and a vertical vector of magnitude Vr. Magnitudes Vr = Vcos a and Vr = 
sin a may be called respectively the horizontal and vertical components of the 
vector V. 
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Figure 2C. An Orbital Transfer. For a satellite moving with velocity V in the 
eiiipticai orbit at position (1), the veiocity is parailel to the oroit at that point and 
can separated into vertical (radial) and horizontal components. At position 
(2), apogee, the vertical component of velocity is zero. If, at that point, an ap- 
propriate horizontal impulse Is provided, the satellite may be transferred to a 
circular orbit with new velocity V' as shown in position (3). In this case the 
magnitude of V' - V = aV Is the change in veiocity. V„ = Vcos r and = 
V sin 7 . 
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Use (32) to find p. 
Use (19) to find e. 
Use (42) to find fp 


Exercise C.3A Calculate the velocities at the points of coincidence for 
each orbit in Exercise B.4.3. The difference in velocities is the Av ra- 
quired for transfer from one orbit to the other when the burns are 
horizontal. 


'.ercise C.3.5. At MECO, the orbiter has a velocity of 25,731.68 ft/sec, 
an aititude of 55 n.mi and a fiight path angie of 1.4° reiative to the 
horizon. 


a) Find Vr and Vr. 

b) Find a, e. 


ANS: 628.68 ft/sec 
25,724.0 ft/sec 

ANS: 3498.96 n.mi 
.0244 


c) Find rg, veiocity at apogee. 

ANS: 3548.3346 n.mi 
25,111.0847 ft/sec 


The orbiter is to be put in a circuiar orbit with a semimajor axis equai to 
rg by burning the QMS engines at apogee. 

d) Find Av required. 

ANS: 312.07 ft/sec 

e) If the orbiter weighs 246,300 lbs, how mucli 
Hydrazine wili be used? 


Exercise C.3.6. Suppose only the initial circular orbit and the final cir- 
cuiar orbit are known. Explain hovi' to calculate the semimajor axis 
and eccentricity of the transfer orbit. 

Hint: The perigee of the transfer orbit must be the same as the radial 
distance of the lower orbit. Similarly, the apogee of the transfer orbit 
must be the same as the radial distance of the upper orbit. 
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C.4. RADIAL BURN TRANSFER 


Transfer from a circular orbit to an elliptical orbit or from one elliptical orbit 
to another may be accomplished by a norizontal burn. The final orbit depends 
on both the magnitude of the fina! velocity and the flight path angle at burn. 

A radial burn from a circular orbit, one directed up or down, generally causes 
a satellite to move to an elliptical orbit. Results of a radial burn for upward 
thrust is Illustrated in Figure 3C. If another circular orbit is desired it could be 
accomplished from the elliptical orbit in the manner discussed earlier. The 
elliptical orbit then would be referred to as the transfer orbit. Note the velocity 
addition loi Figure 3C. Note also the reversibility of the addition. 

Had the satellite been in the elliptical orbit initially with velocity V„ew at posi- 
tion (2), firing the rockets upward to create a downward thrust, or Vr downward, 
would have produced the resultant velocity V, as shown in Figure 4C: The 
resultant find orbit would have been the circular orbit shown. 


Exercise C.4.1. The orbiter is in a circular orbit with an altitude of 160 
n.mi. A radial velocity of 120 ft/sec upward is imparted by firing the 
OMS engines directly toward the center of the earth. 

a) Discuss the change in the orbit. 

Hint: Note that the horizontal component of velocity is not affected by 
a burn at right angles to it. The horizontal velocity is the same as the 
circular orbit velocity at 160 n. mi. altitude (Eq. (71) ), The new total 
velocity is found from (74). The new semimajo'’ axis is found from (16), p 
Ir from (32), the new value of e is found from (19), and o is found 
fron'> 57). 


ANS: 


^circ “ 

25353.94945 

V = 

* new 

26354.2334 

AV = 

0.28398 ft/sec 

®new “ 

Sold -f- 490.50 

®new “■ 

0.004733 

- 

r,„ ^17.1385 

Tp = 

r,|, - 17.056495 


b) How much Hydrazine is burned, if the orbiter weighs 200,000 
lbs? 
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Figure 3C. Radial Burn Transfer. A s iteilite in circular orbit with velocity Is 
shown In position (1). At position (2) engines fire creating a radial thrust giving 
the satellite a radial velocity Vr. Since the velocity V is horizontal and unaf- 
fected, the new velocity is the vector sum of V and Vp with magnitude \/^ 
=^dro+ Ifs direction is the flight path angle 7 , at that point, for the new 
elliptical orbit to which the catellite has been transferred. When the satellite is 
in apogee position (3), an appropriate horizontal thrust could cause it to bo 
transferred from the elliptical orbit to a circular one of radius r^. 


orbit 
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(b) = V„ew + (-Vh). 


Flgure 4C. Vector Addition, (a) parallelogram method, (b) graphic method in 
which vectors to be summed are constructed taii to head with the resuitant 
vector drawn from the tail of the first vector to the head of the last. Note that 

V = V cos 7 , = V sin 7 and 

oirc n«w ' ' R ’new^ " ' i ’'new circ ^ * R' 
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Exercise C.4.2. Suppose the 120 ft/sec increase in speed is imparted by 
a horizontisi burn to the orbiter in a 160 n.mi circuiar orbit, instead of 
the radial burn as In the previous exercise. 

a) Discuss the change in the orbit. 

Hint; Note that the radial component is zero before and after the b jrn. 
The new total velocity is 

= Vciro + 120 ft/sec 

AV = 120 ft/sec 

-f 34.52 n.mi. 

.009444 

rm,o + 68.881 n.mi 
r„|„ - 956.836 ft 

b) How much Hydrazine is burned, if the orbiter 
weighs 200,000 lbs? 

Exercise C.4.3. Discuss the relative merit, in terms of fuel used and 
total change in orbit, of the radial burn and a horizontal burn. 



! Exercise C.4.4. Suppose that the orbiter is in one circular orbit and 

must be transferred to another circular orbit at a new altitude. Is there 
any way to do this with only one short burn of the OMS? 

! 

Exercise CAS. Calculate Ihe radial velocity in an orbit with the 


1 characteristics 

1 

a = 

13,185.3858 

1 

e = 

.7267 

at values of 


0 = 

0”, 180° 

i 

0 = 

90°, 270° 

f 

0 = 

7°, 187° 


0 = 

62°, 242° 




8 ? 


i 



Tabulate the velocities and compare for the cases where the true 
anomaly differs by 180°. 


Exercise C.4.6. Calculate the radial distance, flight path angle, and 
radial velocity for the new orbit resulting from the radial burn in Exer- 
cise C.4.2. 


a) at the point of burn 

b) 180° around the orbit 


Exercise C.4.7. What is the total A v required for placing a satellite In 
the upper circular orbit of Exercise B.4.3. (see Exercise C.3.4.) from a 
160 n.mi circular orbit? The first stage empty case and the second 
stage and the payload weigh 7,436 lbs. How many lbs of a fuel with an 
Ijp = 260 secs is required In the first stage? If the second stage empty 
case and payload weigh 3300 ibs, is there enough fuel in the second 
stage? What is the period of the transfer orbit in hrs, minutes, 
seconds? (Note that ati three orbits are in the same plane. In the real 
problem, we must change planes.) 


Calculation of the time required for an orbiting body to go from one point in 
its orbit to another is important. Position gind time can be accurately related. 
The direct time problem is to find the time associ^teo with a given value of 0 , 
where 0 is the true anomsfy. This is solved as follows. For a given a, e, and 
0 , we may find E by Eq. (58): 


E = 2 tan-' ( ^ - e 


+ e 


tan 


± 

2 


The value of M then may be determined using Eq. (87), 


M = E - e sin E 


Next, we find t/ from Eq (3), 



Finally (t - to) is obtained from the relation 
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where to is the epoch time and t the time when the body reaches the point in its 
orbit defined by 0 . if the epoch is at perigee or apogee, Mo - 0; otherwise, Mo 
must be calcuiated. 


Exercise 0.4.8. Find the time required to travel from perigee to a radial 
distance. 


r = 6800 n.mi, 

for the transfer orbit of Exercise C.4.7. 

Hint; Find 0 from Eq (57). 

ANS: (t - to) = 1^24"' 44* 
where epoch is time of perigee passage 


Exercise C.4.9. Check the answer In Exercise C.4.8. using the calculator 
program of Exercise B.4.1. 


Exercise C.4.10 Find the time required to travel from r = 6800 n.mi to r 
= 18,300 n.mi in the transfer orbit of Lxercise C.4.7. 

Hint: Using Exercise C.4.8., the epoch is 1*' 24'^ 44= and the mean 
anomaly at Epoch is the mean anomaly calculated: Mq = 48” 11’ 40”. 


Exercise C.4.11. The ignition time, direction of burn, and magnitude of 
burn for the second stage is to be determined from radar tracking data 
acquired during the first 1^ 25'^ of travel after first stage ignition. Ap- 
proximately three hours will be required to do calculations. Will there 
be enough time before r = 18,300 n.mi.? 


^erclse C.4.12. A solid rocket, once ignited, cannot be throttled back 
or shut down. Suppose you have to use a two stage solid rocket to 
transfer a payload from the orbi'«er (160 n.mi circular) to another circular 
orbit at 750 n.mi, and the first stage Is much too powerful for the job. 
Do you think you could design a trajectory that would “waste” the right 
amount of fuel, so that delivery can be accomplished? If not, continue. 
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C.S. ENERGY MANAGEMENT: DIRECTIONAL BURNS 

Since fuel economy Is a vital factor In many space misclons, a natural and 
important phase of mission design Is to investigate the conditions that will 
lead to a minimum fuel or energy transfer between the initial and final orbits. A 
transfer of orbit problem rpay be regarded as a change of energy problem. 

First determinations of minimum energy transfer orbits and calculations of 
mission times was done by a man named Walter Hohmann (In 1925!). Let us 
consider two possibilities for transfer between circular orbits. The first, shown 
in Figure 5C, is called a bi-elllptic transfer; three burns are required. The sec- 
ond, a two-stage transfer shown In Figure 6C, has been reviewed in the exer- 
cises and will be referred to as a Hohmann transfer. While the two-burn 
Hohmann transfer will in most cases of interest require less energy, it can be 
shown that the bl-elliptlc transfer Is most energy efficient for cases in vi/hich 
the radius of the final circular orbit is more than sixteen times the radius of the 
initial orbit. But in any case it is much more time consuming. 

Suppose that the arbiter is in a 160 n.mi circular orbit, and an upper stage is 
deployed for transfer into a geocentric orbit. The transfer is minimum energy 
(Hohmann transfer) with the characteristics: 

Basic 160 n.mi orbit; 
a = 3603.930885 n.mi 

V = 25353.9497 fps 

P = 1^’30"’27* 

Transfer orbit: 
a = 13,185.3858 

e = ,7266723 

Vp = 33,315.8263 fps 

V« = 5,273.8074 fps 

P = lO'' 32"' 56" 

Geosynchronous orbit: 
a = 22,766.84089 n.mi 

V = 10,087.4748 fps 

P = 23" m"' 04' 
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(3) 


final orbit 





Figure 5C. Bi-elliptic Transfer. The first of three horizontal burns is made at (2) 
resulting in velocity V,. The second burn is at apogee (4) of the elliptical orbit, 
resulting in velocity V,. The third burn, at perigee (6) in the new eliiptical orbit, 
is made to give an opposing velocity in order to effect the final circular orbit. 




Figure 6C, Hohmann Transfer. This is a minimum energy transfer between two 
circular orbits, accomplished by two horizontal burns. The first burn occurs at 
(2) imparting velocity AV,; the second burn occuid at (4), 180 ° away, imparting 
velocity AVj. Vj - V = AV, + aV, = AV. 
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You may verify the parameters for the transfer orbit using previous exercIseSi 
and that the required AV for the two upper stages is given by 

1st stage AV - 7961.68 fps 

2nd stage 6 = 4813.67 fps 

Suppose now that the upper stage has two motors, a solid rocket for the first 
stage and a liquid fuel rocket for the second stage. The solid rocket cannot be 
throttled, and has a fixed Av of 8,200 fps while the liquid stage can be throttled 
to any desired AV. The prescribed Hohmann transfer orbit cannot be made 
since the AV is too large. But now consider the problem of choosing a direc- 
tion of first stage burn such that excess Av in the first stage is wasted, and of 
choosing a direction and value of the second stage burn to accomplish the 
desired circularization at geosynchronous altitude, as illustrated in Figure 7C. 

The procedure for solution Is as follows. Choose the burn direction for the 
first stage such that the horizontal component of velocity gained is the Av re- 
quired to go from the circular orbit to the prescribed transfer orbit, as shown In 
Figure 8C. Using the orbital characteristics add given information, we see that 
In two triangles of Figure 8C, 


Vc = 25353.9497 ft/sec 


AV^ = 79618766 ft/sec 
Vciro + = 33,315,8263 ft/sec = Vh, 


1st stage AV = 8,200 ft/sec. 


The angle P is given by 


P - COS gjgyg 


P = = 13.8417 » 


From Figure 8C, 


Vr = (1st stage AV ) sin p . 


Vr = 8200 sin (13.8417°), 



Figure 7C. Orbital Transfer Away From Apogee. At (1) a satellite in circular or- 
bit with velocity Vq. At (2) engines are fired producing a velocity AV,. The new 
velocity Vt is the vector sum of AV, and Vq. At (3) engines are fired in a 
predetermined direction producing a velocity AVj. The resultant velocity of 
the satellite will be the vector sum of aV, .and Vj at that point. If the sum is 
equal to Vq^o and directed horizontally, the satellite will be placed in the 
geosynchronous orbit shown. 
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Figure 8C. Velocity Addition for Directed Burn. First stage Av Is vector 
sum of avreo and Vp,. Angle p is the burn angle, The resultant velocity after 
burn, V„^, is the vector sum of the Initital velocity In circular orbit, Vc, and the 
first stage Av, y Is the flight path angle for the transfer orbit. The horizontal 
velocity after burn, V,,, Is the sum of and aVreq, 


V„ - 19eiJ648 ft/sec 


As a check try 


=.y(1st stage AV)»-( AVpeq)= 


Vfl = V(8200)* - (7961 .877)^, 


Vr = 1961,7648 ft/sec 


The new orbital velocity is 


= vv?- + n 


V,»w = V(33,315.8263)2 + (1961.7648)2, 


V„^ = 33,373.53447 ft/sec. 
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The flight path angle for the first stage is found from Eq, (65). 

7 = tan”' ^ 


7 


tan*' ( 


1961.7646 

35,315.8263 


7 = 3.36990 *, 

The new transfer orbit has a set of characteristics different from that given. 
The parameter p Is obtained from Eq. (32). 



3603.93088 33,315.8263 

^ ^ 62750.278V ' 6076.1155 

p = 6222.8079 n.ml. 

The new semimajor axis is given by Eq. (16). 

a L 

2-rJ^ 

a = 13480,6653 n.ml. 


The new eccentricity is given by Eq. (19). 


e 



e = 0,73375074 

The line of apsides of the new transfer orbit has been rotated by the radial 
burn. The new true anomaly is given by Eq. (54). 
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© = 


tair^ ( 


p tan '' 
P- r 


) 


© = 7,9649“ 

We have now completed the calculation of 
(I) the required burn direction, 

(ii) the characteristics of the new transfer orbit. 

Th») direction and magnitude of the 2nd stage burn is determined by the in- 
tersection of the transfer orbit and the required geosynchronous orbit. To 
determine where this intersection occurs (see Figure 7C) we determine the true 
anomaiy for r equal to the radius of the geosynchronous orbit. Using Eq. (57). 


e 


cos-’ ( 


P - ^OEO 

6fGE0 


„ . , , 6222,8079 - 22,766.8409 , 

6 = cos-' ( 

' (0,73373) 22,766.8409 ' 


e' = 172.0355“. 

Note that the angle of intersection, e", when added to the angle at burn, is 
180“. Under the circumstances of burn we have chosen, orbital intersection will 
occur at angles d' = 180“ - 6 and at d' = 180“ + e . 

Assume that we have chosen for orbital change the second intersection 
point, e' = 180“ + e , The radial velocity at this point of intersection is given 
by Eq. (77). 


Vr 


sin (180 -f- 0 ), 
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where 0 le the true anomaly of the 1st burn. Inaerting the data, 

= 1961.7661 ft/eec. 

Note that this differs from the radial velocity of the first burn by: (I) It Is In op- 
posite direction; (II) It differs In magnitude by less than 0,002 ft/sac, which we 
can safely attribute to rounding error In the computations, See Exercise C.4,5, 
The orbital velocity Is given by Eq, (70), 

V ^ (t4) 


or, 


V = 5,626,861 ft/sec. 


The flight path angle Is given by Eq, (63), 

7 = sin-’ 

V 


yielding 


7 = 20,4043”, 

The horlzofjtal velocity Is given by 


= V cos 7 , 


or 


V„ = 5,273.307 ft/sec 


As a check on the calculated horizontal velocity, use 


to obtain 


V„ = 5,273.8068 ft/aec, 

Knowing the velocity required to circularize at geosynchronous altitude 
to be 


Vqeo = 10,087.4748 ft/sec, 


The aV requlr, * Is 


Vqeo • Vh = 4,813.668 ft/sec. 

The second stage must be burned in a direction with a magnitude that (i) 
cancels the radial component In the transfer orbit, and (ii) provides the A v 
necessary in the horizontal direction. From Figure 9C, the magnitude of the 
2nd stage Av required is 


2nd Stage Av = 

= ^(481 3.688)* + (1961.7661^ 


2nd Stage av = 5,198.069 ft/sec 
The angle at which the 2nd stage must be burned is given by 


(3 = 180 - (90 -)t|) - cos-' ( 




2nd stage A v 


) 




= 90 + l'y| - cos-' ( 


1961.7661 
'5198.069 ' 


where we have used the facts that the sum of the angles in a triangle is 180 
and that the total on a straight line is 180°. Then 



Figure 9C. Second Stage Velocity Addition, The resultant of V^,^ and the 2nd 
stage av Is the final velocity for the desired circular orbit. This vehclty 
I must be horizontal and equal to Vqeo- The radial component of Av must be 

identital to and cancel the radial component of V^^. The flight path angle y Is 
Identical to the flight path angle caused by the first burn, 
i 
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(3 = 42.5772 


and the burn direction has been found, completing the calculations. Note 
again that In the real world, a change of plane, which we have Ignored, is also 
necessary. 


Exercise C.5.1. Since an intersection with the transfer orbit Just deter- 
mined occurs at a true anomaly of 


100-0 


in the new transfer orbit, derive the magnitude and direction of the 2nd 
stage burn required to circularize into a geosynchronous orbit at this 
first Intersection. 0 is the true anomaly at 1st stage burn ( 0 - 
7.9649“). 


Exercise C.5.2. Find the mean anomaly at epoch and the epoch after 
time of perigee passage for the new transfer orbit of the last exercise 
at © = 7.9649“, 


t„ = o'' 01"' 31“ 

Mo = 0.8333“ 

Verify these values with the calculator program of Exercise 8.4.1. 


Exercise C.5.3. Find the times required to Intersect the geosyn- 
chronous orbit altitude In the new transfer orbit, after the time of first 
burn. 

ANS: 

1st opportunity M, = 145“ 22’ 19" 

(t - U = 4" 21'" 11“ 

2nd opportunity Mj = 214“ 32’ 31" 

(t - to)j = 6'’ 28"' 34“ 

These times are lapsed times after the first stage burn. 
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Exercise C.5.4. Calculate the Av required for escape from earth for a 
160 n.mi altitude orbit (Eq. (79) ), Compare this with the total av re< 
quired to place a payload in a geosynchronous orbit from 160 n.mi, 
Compare the two. 


Exercise C.S.S, The orbiter is in an orbit with the characteristics 

perigee altitude = 120 n.mi 
apogee altitude = 260 n.mi 

An upper stage is to be deployed from perigee for transfer into geosyn- 
chronous orbit. The first stage is a solid rocket with a av of 8,200 
ft/sec. The second stage is a liquid fuel stage. 

Design the upper stage burn directions, magnitudes, and burn loca- 
tions. There are two opportunities for upper stage ignition during each 
revolution In the transfer oibit Find the times of each opportunity 


C.6. EVALUATION 
C.6.1. Learning Objectives 

The principal objectives of Section C were fourfold: first, to show how sets of 
parameters characterize an orbit and to develop methods for calculating the 
position and velocity of a satellite in the plane of its orbit; second, to show how 
particular orbits may be attained through orbital transfers; thirdly, to show how 
orbital maneuvers may be influenced by upper stage rocket characteristics and 
fuel efficiency requirements; and lastly, to raise the level of analysis from 
single to multistage and multivariable. 

Having completed your study of this section and related references, you 
should have an understanding of 

• the elements required for describing orbit characteristics 
and satellite position in the plane 

• the problems relating to transfer orbits and fuel efficiency 

• the use of burn directions, thrust magnitudes, and orbital 
position in effecting orbital transfers 


I 
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C.e^ Outeomo Measures 


You should now be able to 


1. Use the appropriate basic transfer equations in solving prob- 
lems, 

2 . Use the appropriate formulas and srecific date to calculate 
(a) the position of an orbiter, (b) the velocity of an orbiter at a 
particular point In Its orbit, and (c) the time required for an or- 
jjiter to travel frrm one point to another. 

3. Describe orbital transfers. 

4. Add and resolve vectors Information components; 

5. Design upper stage burn directions, magnitudes, and burn 
locations for designated orbitai transfers. 

6„ Define fuel •■equlrements for orbital maneuvers, given appro- 
priate information. 


D. TIME AND REUTIVE ORBITAL POSITION 


D.O. INTRODUCTION 

We have discussed the position and velocity of a satellite, arbiter, or other 
spacecraft relative to the apogee and perigee ooints of Its orbit, and orbital 
transfers within the orbital plane; now we turn our attention to a more com- 
plete description of the position of an orbiting object, one relative to Its base — 
the earth — and its energy source — the sun. The earth Is a ball, not quite 
spherical, traveling at high speed In Its orbit about the sun and spinning like a 
top about Its axis. When we calculate the velocity of a satellite, It is usually the 
velocity relative to an earth <;/urface, and we want to be able to specify its posi- 
tion relative to some position on that surface. 


0.1. CELESTIAL COORDINATES 

The position of stars and other celestial or space objects may be specified 
in the same way that a position is defined on the earth’s surface. The location 
of any spot on earth Is given by latitudinal and longitudirial meridians. Latitude 
is measured north or south of the equator, which is assigned a value of 0“, as 
shown in Figure ID. The po/es are to the north and south, Longitude \s 
measured east or west of the Greenwich Meridian, which is the 0° reference. 
One nautical mile at the equator was originally the equivalent of one minute of 
arc. 

To define a set of Celestial coordinates, we construct a Celestial Sphere 
with the earth at center as shown in Figure 2D. The Intersection of the pro- 
jected plane of the earth’s equator with the Celestial Sphere defines he 
Celestial Equate' A line from the Celestial North Pole to the Celestial South 
Pole passing th.ough the point in the sky, the First Point of Arles or Vernal 
Equinox, where the sun crosses the equator in the spring, provides the 
longitude type reference for the sky. The symbol, j , sign of the Ram, is used 
to denote the Vernal Equinox and may be recognized by those of you born 
under it. Angles are measured to the east along the Celestial Equator from the 
First Point of Aries. 


Exercise D.1.1. Using the value of 1 n.nii = 1 min of arc, calculate the 
equatorial radius ir< nautical miles. Compare this with the value in Ap- 
pendix A. (See, for example, me entry nautical mile in Webster’ » 
Unabridged Dic tionary for an indication of the reason tor the discreparr 
cy.) 
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Celestial North Pole 



Figure 2D. Celestial Coordinates. The line from the Celestial North Pole 
through the First Point of Aries to the South Pole provides a Celestial 
"longitude” reference. An object positioned as shown is located through the 
angle RA, called the right ascension (defining "longitude”) and the angle I, the 
declination (defining "latitude”). 


I 


Longihide 


Greenwich, 


rarnr 


7A 


.^. 30 ° 


\ 


Greenwich 
Meridian (0°) 


Figure ID. Earth Meridians, The Greenwich meridian at 0° is longitudinal 
reference for the earth, The equator is the latitudinal reference at 0°. The 
angle shown in the equatorial plane Is called the Greenwich Hour Angle (GHA). 
GHA = GHA^ + oj * (t) Is the rotation rate of the earth and t is the el lapsed 
time from epoch. 



0^ TIME MEASUREMENT 


The clock time Interval between euccettslvti i warl Jian passages of the sun Is 
called a truo solar day. That is to say, If a time check is made when the sun Is 
directly ovedread on one day and again when the sun Is directly overhead the 
following day, that time interval is the true solar day, as illustrated in Figure 3D, 
The length of the true solar day varies throughout the year, since it is depen- 
dent upon the speed of the earth which varies (according to Kepler's second 
law) In Its orbit about the sun. The average cI(k:K time length of the solar day 
over the year is called the mean solar day. The mean solar day Is divided Into 
24 hours or 86,400 mean solar seconds. 

We know the motion to be relative. It is accepted fact that the earth spins on 
Its axis and moves In orbit about a focal p'jint at the center of the sun. But 
such motion can be viewed from differcynt perspectives (different coordinate or 
reference systems). Choosing to regard the earth as stationary in its orbit we 
obtain a different view: the sun would appear to move in orbit about the earth 
with a period of one year. Choosing to regard the earth as stationary in Its rota- 
tional motion, the sun would a,ppear to orbit about the earth with a period of 
one day: this Is oor normal ‘‘illusion, ' with the solar day defined In terms of 
successive earth n>')r!dlan passages of the sun In Its orbit, shown In Figure 4D. 
We use the view or system that Is most convenient for our problem description. 

The time divisions of the clock for a mean solar day define the quantity 
called Universal Time (UT), They are also the same for Greenwich Mean Time 
(GMT), 12 hours GMT (noon) defined to be the Instant of passage of the mean 
sun over the Greenwich Meridian (0® longitude). The Local Mean Time (LMT) at 
any longitude on earth can be calculated by 


or 


LMT = GMT 4- (East longitude)/ 15 
LMT = GMT - (W'ast longitude)/15 


where the longitude is given in degrees and fractions of degrees. 

The time Intenral between successive meridian passages of the vernal 
equinox is called a sidereal day. It is roughly the time for one earth rotation, as 
was shown in Figure 3D. But because of the “wobble" of the earth's poles, 
called precession, the sidereal day varies slightly in length. As with the solar 
day, we may define a mean sidereal day, Greenwich Mean Sidereal Time (GST), 
and Local Sidereal Time (LST). A sidereal time interval is somewhat shorter 
than a solar time interval: 

sidereal time interval = 0,997269^84 
(solar time interval) 
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Figure 3D. The Solar and Sidereal Day. With the earth considered in its mo- 
tions, the time Tg at P is noon with the sun directly overhead in (a). When the 
earth makes one complete rotation, the point P will be located as shown in (b). 
Due to the movement of the earth in its orbit, the sun is now directly overhead 
at the point P’. The time required for the earth to rotate from P (a) to P (b) Is one 
sidereal day. The time required for the earth to rotate from P (a) to P (b) is one 
solar day. If P Is on the Greenwich Meridian, {the time at P (a) Is 12*' Greenwich 
Mean Time (GMT). 
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ExerelM D.2.1. Find the length of a sidereal day in terms of the mean 
solar day. 


Exeteise D.2.2. A point near Houston, Texas, has the longitude 95” 00’ 
28"W. Find the Local Mean Time at 12^ Greenwich Mean Time, 


Exercise D^.3. Find the length of a mean solar day in sidereal time. 


D.3. ANQUUR TIME MEASURE 

The time, in angle measure, since the last passage of the first point of Aries 
over the Greenwich Meridian is called the Greenwich Hour Angle (GHA) (See 
Figure 1D). The GHA expressed in time measure is identical to the Greenwich 
Mean Sidereal Time (GST). Thus, conversion of time measure to angle 
measure is based on exactly 15 seconds of arc per sidereal set;>nd, or 

15.041069 sec of arc = 1 mean solar second 


Exercise D.3.1. Using the number of revolutior>s per day from Exercise 
A.7.I., calculate the number of seconds of arc per mean solar day, and 
compare your answer with the number above. Calculate the number of 
degrees per solar hour. 


The GHA of Aries can be computed at any time from the equation 
GHA (t) = GHA (to) + (t - to) 

where time is in mean solar seconds and GHA (t) is in degrees and fractions of 
degrees. 


Exercise 0,3.2. Compute the GHA of Aries at 17^' GMT on July 2, 1979, 
using the epoch value, GHA (to) = 278° 27’ 12”, on July 1, 1979. Com- 
pare your answer with the A.E Nautical Almanac entry. 


D.4. TIME CORRECTIONS 

Tims is given in the American Ephemerls and Nautical Almanac in what is 
called Ephemerls Time (ET). This is a uniform time measurement based in prin- 
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cipal on observations of the planets and moon. There is no theoretical connec- 
tion between Ephemerls time and Universal time, or GMT, which fluctuates as a 
result of changing earth motion. The difference is observed over several years 
and projected forward one year at a time. For example, the predicted difference 
between GMT and Ephemerls time In 1979 Is given as; 

Jan. 1, 1979 -r-49.9 seconds 

Apr. 1, 1979 -*-50.1 seconds 

July 1, 1979 -H50.4 seconds 

Oct. 1, 1979 -r-50.7 seconds 

The difference is quite small, but correction to Universal Time should be made If 
high accuracy is required. 

For example, the right ascension of the sun at 0^ July 1, 1979, is given in the 
Almanac as 6*' 37"' 23.45*. The correction for GMT Is 

X X 15.0 = 2.16" 

86400 

Therefore, the corrected right ascension for o'* GMT on the given date is 99° 20’ 
54.16". 


Exercise D.4.1. Use the American Ephemerls and Nautical Almanac to 
find the Universal Time of the meridian passage of the equinox at 
Greenwich on July 1, 1979. 


0.5. THE RIGHT ASCENSION AND LONGITUDE 


The angle measured to the east along the Celestial Equator from the first 
point of Aries to the celestial meridian of an object in the sky is called the right 
ascension (RA) of the object. This has been Illustrated in Figure 2D. It is the 
difference between the Greenwich Hour Angle of Aries and the GHA of the ob- 
ject. 


The Local Hour Angle (LHA) is the time since the local meridian passage of 
the first point of Aries. It is identical to the Local Sidereal Time when ex- 
pressed in time measure. 


LHA(T) = GHA(t) 


-f East longitude ( X E) 
- West Longitude ( X W) 
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If the right ascension of a celestial object Is equal to the Local Hour Angle, i.e., 

when RA = LHA(T) 

the object is directly over the local meridian. Therefore, the earth fixed 
longitude, \ , of a subsatellite point can be found if the right ascension of the 
satellite can be found. The earth iongitude is given by 

= RA-GHA(t) 

where measurement is in the eastern direction (see Figure 5D). 


Exercise D.S.I. Using the Nauticai Aimanac, find the right ascension of 
the sun at 0*' GMT on Juiy 1, 1981. 


Exercise D.5.St. At 0^ GMT Juiy 1, 1981, caiculate where on earth it is 
exactiy noon. Note that it is the “subsateiiite” iongitude of the sun re- 
quested. Compare your answer with the GHA (sun) from the Aimanac 
for that date. 


D.6. THE RiQHT ASCENSION OF THE MEAN SUN 


in our work a rapid means of caiculating the right ascension of the mean 
sun, RA (sun), is an advantage over the use of tabies, for reasons of con- 
venience and efficiency. The orbitai eiements of the sun are used; again, we 
pretend that the sun is in orbit about the earth when doing these caicuiations. 
The coordinate system used is the true equator and true equinox of date, as 
shown in Figure 6D. 

Consider the epoch 0*’ GMT, January 1, 1979 (JD 2443874.5). The orbital 
eiements of the sun will be 

g = 282.5711432” -^ (0.0000470684”) d, 

M = 357.47863 (0.985600267) d, 

e = 0.016736825, 

G = 23.44207294". 

g is called the argument of perigee, the angle measured in the 
plane of the ecliptic (sun’s orbital plane illustrated in Figure 
8D) from the true equinox of date (First Point of Aries) to the 
perigee of the sun'.s orbit, in degrees. 
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Figure 5D. Longitude of Subsateilite' Point. (GHA (t ) is the angie measured 
from the Greenwich Meridian to a iine pointing to the First Point of Aries or the 
Vernai Equinox. RA is the right ascension of the sateiiite, measured in the 
equatoriai plane from the Vernal Equinox as shown. The satellite Is directly 
over the point P which lies on the earth longitude X . With GHA (r) greater 
than RA, the longitude Is west with X„ = GHA (t) - RA. 
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North Pole 



Figure 61). The Ecliptic. The yearly path that the sun appears to follow among 
the stars obsen/ed from earth is called the eliptic. It makes an angle e =23° 
27 "' with the celestial equator. Point A is the Vernal Equinox or First point of 
Aries, the ascending node; with the sun in that position, on March 21, the 
length of day and night are equal. Point C, the decending node, occurs on 
September 21. The coordinate system is the True Equator and Equinox of 
Epoch (TEE) reference axis system. The epoch Is defined to be 0^ 0"' 0“ GMT on 
the user specified base date. 
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Figure 7D. Location of Orbiting Body. A body in orbit is specified by its right 
ij ascension, RA, and deciination, measured in and from the equatorial plane. 

I The angle u is measured in the orbital plane. 

I 
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equinox 



Figure 80. Angles in the Plane of Motion. P refers to the orbit perigee and A to 
I apogee. The angle g is the argument of perigee, 6 is the true anomaly, and U, 

i measured from the ascending node to the object, is the argument of latitude. 


I 

I 


t 
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d is the number of Julian days elapsed after 0** GMT, January 
1, 1979, In days and decimal fractions of days. (Note; The 
correction term added to g represents the perturbation of the 
earth’s orbit by the moon and planets. It amounts to 0.1694 
secs of arc per day rotation of the line of apsides of the 


M is the mean anomaly of the sun In its orbit. The term involv- 
ing the number of Juiian days since epoch reflects the 
period of the sun In Its orbit; 


P = 


360 deg 
0.985600267 


deg/deg = 365.259641 days 


e is the eccentricity. It chan;/es due to perturbations of the 
moon and planets about 41,8 parts per million per century 
(decreasing during our geological era). We will ignore this 
change. 

6 is the obliquity of the ecliptic (the inclination of the plane of 
the sun's orbit) in degrees. The obliquity of the ecliptic is 
also changing at the rate of 47 secs of arc per century. We 
will Ignore this change. 


The procedure for calculating the right ascension is as follow.s. 

(a) Find the Julian date of the new date. 

(b) Subtract the Julian date of epoch o'' GMT January 1, 1979 
(JD 2443874.5) to find d, 

d = JD - 2443874.5 

Note that the Julian date may include a decimal fraction of a 
day. For example, at 17^ 32™ 03.3* GMT on 30 June 1979, the 
Julian date is 

2444054.5 + = 2444055.23 

24 

(c) Find the new elements of the orbit at the new date. 

(d) Find the true anomaly (see programs written in Exercises 
B.2.1 and B.4.I.), 

0 = M -f 2 e sin M -r- -f- e* sin 2 M 

4 

-I- ei (13 sin 3M - 3 sin M) 

12 
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(e) Find the argument of latitude, u, defined as the angle In the 
plane of motion from the ascending node, i.e., the true 
equinox of date to the position of the object), 

u = g 0 


(f) Calculate 


X = cos u, 
y = cos e sin u, 

RA jSun) = tan-' (-—) 

(g) Do a quadrant check. 


Exercise D.6.1. Repeat the calculation for the Julian date for 17*' 32”' 
03.3* GMT on 30 June 1979. Rearrange the computations as necessary 
to carry the number of decimal fractions of a day to 0.1 second of time. 
How many decimal places must be retained? 


Exercise D.6.2. Find the right ascension of the sun at 17” 32'" 03.3’ 
GMT on June 30, 1979. In your solution you should find: 

d = 180.73059 days, 

g = 282.57965“, 

M = 175.60675, 

u = 98.18640 + 57.29578 (0.0025118), 

W = 98.18640 + 0.143911, 

X = -0.144880, 

y = 0.907783, 

RA(sun) = -80.93222“ + 180“, 

RA(sun) =99“ 04’ 04”. 

Exercise D.6.3. Check the result of the last exercise against the tabular 
entry in the American Ephemeris and Nautical Almanac. What is the , 
discrepancy? 
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ANS: at O^ET, RA = 6** 33"’ 14,94’, correapondlng to 90.3125* 

at + I?” 32"’ 03.3* ET, RA = 99* 04' 08" 

(0.130590 X 248.51 = 181,5589 = 2723.385 sac arc = 
0.756496*) 

ET to GMT correction, RA = 99 * 04’ 10.16" 

<501 X X 15 = 2.16 sec of arc) 


Exercise D.6.4. (a) What clock error corresponds to 6.2 secs of arc error 
In RA (Sun)? 


ANS: 143.7 sec 

(b) How many decimal places In the Julian Day number does this repre- 
sent? 


ANS: IE -06 


Exercise D.6.5. What Is the maximum error, 6 , in RA (Sun) for an 
observer on the equator If he uses the tabular entry in the Alm.anac? 
See Figure 9D. 


ANS: 9 secs of arc 


Exercise D.6.6. What is the maximum error in RA (Sun) for the orbiter in 
a 160 n. mi orbit using the tabular entries In the Almanac? 


Exercise D.6.7. For an obsen/er at earth’s distance, what is the max- 
imum error in RA (Sun) If one limb is observed (see o In Figure 9D) in- 
stead of the center of the sun? Note that th,e sun has a diameter of 
752,827.06 n.ml, and find its distance on June 30, 1979 from the 
Almanac. 


ANS: 15’ 45" arc 


Exercise D.6.8. Check this last answer against the apparent semi- 
diameter of the sun given in the Almanac for 0*’ ET June 30, 1979. 
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0.6.8. Write a calculator program for calculating tfie RA (Sun) 
and GHA (Sun) given t!ie time, GMT, and number of Julian days 
elapsed sinch epoch. Choose the epoch to be o'' GMT January 1, 1876. 


D.7. THE DECLINATION OF THE SUN 

The declination of a celestial object Is the angle from the celestial equator to 
the obje^ t, measured positive to the not n, negative to the south. Check 
Figures 2 and 8D. The declination of the sun is given by 

e = sin** (sin e sin u), 

where 

e is the declination, 

C is the obliquity of the ecliptic, 

u is the argument of latitude, calculated as In the previous 
section. 


Exercise 0.7.1. Find the declination of the sun at 17*' 32"' 03.3* GMT on 
June 30, 1979. Compare your answer with the Almanac tabulated 
value. 

ANS; 23“ 10’ 49’»; 23“ 10’ 36” tabulated 


Exercise 0.7.2. Add the computation of the declination to the 
calculator program of Exercise 0.6.0. 


DJ. SUNRISE AND SUNSET 

Daylight and darkness onset times are calculated in terms of the maximum 
central angle of daylight. This was shown in Figure 90. P would be 90“ if the 
sun were a point (infinitely far away), if the earth were a perfect sphere, and if 
the atmosphere did not bend the rays of the sun. The correction for the fact 
that the sun has an obsen/able diameter gives 

= 90“ 15’ 45" 

The longitudes of sunrise and sunset are found from 
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Mhtrs 

! 

LAT 


X «yn^!»o s qha (sun) ± cos-‘ ( ^ ^ ' ) 

sunset ' ^ cos 1 cos LAT 


Is the declination of the sun at date, 

Is the given latitude, 

Is the central angle ( ^ = 90“ 15’ 46"). 


QHA (sun) = GHA (T) - RA (sun) 


gives the longitude on the earth where it Is exactly noon. 


Exwelse DAI. Find the longitude of sunset at 17^ 32"' 03* GMT June 
30, 1979 at latitude 40 “N. 


Let us look at the solution for the above exercise. From the Nautical 
Almanac, the QHA (T ) at 0*’ GMT June 30, 1979 Is 277* 28.0', At 17" 32"* 0,3.*.3, 

GHA(T) = 277, 46667“ + 15,041069(17.53425), 

GHA(T) = 181“. 200534, 

From Exercise D.6,2, RA (sun) = 99 * 04' 04’’. Therefore, 

GHA (sun) = - 82 * 07"’ 58*. 

From Exercise D.7.1, the declination Is 

I = 23* 10’ 49", 


Using the given relation, 

X sunset = -82 *07’ 58" + cos-’ ^-^^^5^^:25301^ 

= -82* 07’ 58" + 111 *27’ 25” 

= 29“ 19’ 27" 

Exercise D.8A (a) Find the time of sunset for the previous exercise in 
GMT. 

ANS: 17" 32'" 03.3* 

(b) Find time of sunset in local time. 
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We will note that the solution for (b) above Is 

12*' + = ig** 25'" 50‘ 

10 

That solution follows because 

(a) the longitude of high noon corresponds to 12^ and 

(b) the longitude to the east converts to time at the rate of 15* 
per hour, later. 

Check the local time against the entry In the Nautical Almanac for the June 30 
date. 

When the second term in the equation for longitude of sunrise/sunset is 
160* or undefined, the sun never rises (or sots) at these latitudes. Thus, the 
condl ‘on for a continuous daylight lor darkness) is 

cos 0 - sin I sin LAT ^ .| 
cos I cos LAT 


This occurs when the absolute value of the latitude Is greater than ^ minus 
the absolute value of the declination, 


UT 



= 


I 


Thus, for a declination of 23“ 10' 49" (17" 32'" 03.3' GMT on June 30, 1979) con- 
tinual daylight occurs above latitude 67° 04’ 56”N and continual darkness 
exists below 67° 04’ 56’’S. Those latitudes must be modified to account for the 
oblate figure of the earth. Check this result (67° 04’ 56") against the 
Sunrise/Sunset table in the Nautical Almanac for this date. 


D.9. OAYLIQHT/DARKNESS FOR A SPACECRAFT 

For a spacecraft above the earth, the central angle ^ must be modified as 
shown In Figure 10D. The central angle is given by 

= COS-’ (^)+ 90° 15’ 45" 

Exercise D.9.1. Find the central angle for the orbiter in a 160 n.ml cir- 
cular orbit. 

ANS: 107° 23’ 57" 


118 






to center of sun 


orbit of spacecraft 


;darkness/ 


Figure 10D. (a) shows a view perpendicular to orbital plane. is the central 
angle for the spacecraft at a radial distance r„ with (3o equal to 90° 15’ 45". (b) 
shows a view parallel to equatorial plane with orbital plane given parallel. 1, is 
the declination of the orbiting spacecraft and 1,„n the declination of the sun. In 
(a) is shown possible sunrise and sunset points, but (b) shows that the 
spacecraft is always in sunlight due to the declinations of sun and orbital 
planes. Latitudes above LAT, on earth can be seen to be always in sunlight. 


1 


The coordinates of the spacecraft and the sun are obtained by projections 
against the celestial sphere, so that the equation for sunrise/sunset becomes 

RA (8) = RA (sun) ± cos-' ( ^ '» W ) 

cos I. cos l.,„ 

where 

1^ is the declination of the sun, In degrees, 

I, is the declination of the spacecraft, in degrees. 

For the special ciise of a spacecraft in an equatorlsl orbit, I, = 0”, the relation 
reduces to 


RA (s) = RA (sun) ^ cos-' ( cos P ) 

I Sun 


Exercise D.9.2. For a geosynchronous satellite in a circular equatorial 
orbit, find the amount of time It spends in darkness on July 1, 1979. 

ANS; 


0 = 171” 33' 43” 

RA (Sun) = 99“ 20’ 54" 
l.„„ = 23” 09’ 40” 

Inverse cosine Undefined- 
Spacecraft always in sunlight 


Exercise D.9.3. Find the amount of time the geosynchronous satellite 
spends In darkness on September 23, 1979. 


For a solution to Exercise D.9.3, we use the program of Exercises D.6.9, and 
D.7.2. On 0*' GMT, September 23, 1979, RA (sun) is found to be given by 


Then 

or 


RA(Sun) = 179” 25’ 40” 

RA(Sun)= 40” 15’ 01” 

RA(s) = 179” 25’ 40” ± 171 ”33’ 56”, 

sunset 

RA (s) = 350 ° 59’ 36” Sunset 
RA(s) = 7” 51’ 44” Sunrise 


It follows that 


angle in Darkness = 16 * 52* 08** , 

total period in hours = 23.93446 

time in darkness = x 23.93446374, 

3du0 

t =; 1*’07"’ 17* 

The ratio of time spent can be found this way only because It is a circular orbit 
and the motion is uniform throughout the orbit (I.e., constant velocity). 


D.10. EARTH FIXED LONGITUDES 

The earth fixed longitude is calculated by the equation given in D.5 (see 
Figure 50). 


X E = RA(object)-GHA(T) 

SImil&rly, If the longitude of an observation and the time of an observation Is 
known, 


RA (object) = X E + GHA(t) 


Exercise D.10.1. A satellite passes directly overhead at Houston (95° 28* 
W) at 14*’ 23'" 16* GMT on July 2, 1979. Find the right ascension of the 
satellite. 


ANS: 40° 22’ 45** 


Exercise D.10.2. The arbiter is to deploy an upper stage in full sunlight 
under the condition that the angle between the sun and the flight path 
angle of the upper stage shall be at least 35° and no greater than 80°. 
The deployment is also to be over a tracking site at 91 ° 32* W. Launch Is 
to be between 10" GMT and 11^ GMT on July 15, 1981. How long after 
the launch will deployment be possible? We will assist in the solution 
starting as you should, with a sketch. 
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The RA (Sun) must be between the values 
RA (Sun) = 91“ 32' 0” + 10° + GHA (T) 
and 

RA (Sun) = 91° 32’ 0" + 55° + GHA (t) 
The required GHA (T) at deploy Is between 
GHA (T) = RA (Sun) + 91 ° 32' 0” - 10° 
and 

GHA (T) = RA (Sun) + 91 ° 32' 0" - 55° 
Use Programs of Exercise D.6.9. 

» Find RA (Sun) at lo'’ GMT on July 15, 1981 
• Find GHA (T) at 10" GMT on July 15, 1981 
GHA (T) at 11" GMT on July 15, 1981 
(Note RA (sun) does not change rapidly.) 
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• Find differences between QHA (T) at 10*’ and 11*’ and re- 
quired QHA (n ) 

• Convert to lime 

• Correct for change In RA (sun) at new times 

A'NS: • The deploy window opens at 14*’ 32"’ 0* GMT 

• Deploy window closes at 17*' 32"’ 01* 

• For lO" GMT launch, the times prior to deploy are 4*’ 32"' 0* 
and 7" 32"’ 01* 

• For 11*’ GMT launcli the times are reduced by 1 hour. 

• The above Includes a correction of 0*’ O'" 45* window opening. 

• o" l"’ 16* window closing for change In RA (sun) between 
above times and 10** GMT. 

• Does not Include a correction for the seml-dlameter of the 
sun. 


Exercise D.10.3. For an orbiter In a 1^ n.ml circular orbit, assuming no 
latitude requirements for deploy (e.g., no equatorial crossing display re- 
quirement), how many opportunities, at most, will the orbiter have to 
deploy the payload In the above window? 

ANS: 2 opportunities 

Note: If an equatorial crossing deploy requirement exists (e.g., plane 
change requirement), no opportunity may exist in this window for this 
launch time. 


Exercise D.10.4. In the above window, how much time will the crew 
have to deploy the payload? 

ANS:0*’iri8* 


t 
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Exeieise D.10.5. Suppose something precludes the deployment. Where 
Is the next window? 

ANS: window opens on July 16, 1981 14^ 32"’ 05* GMT 
closes 17*’ 32"’ 04* 

Exercise D.10.6. Suppose the satellite In the previous example Is to be 
placed In geosynchronous orbit by a direct Hohmann transfei. What Is 
the Earth-fixed longitude of placement? (Assume It Is deployed from a 
160 n.ml circular orbit.) 

ANS;9°08'E 

Exercise D.10.7. If a geosynchronous satellite Is to remain over the 
Houston, Texas, meridian (95° 28' W), and It has been placed In geosyn- 
chronous orbit by a Hohmann transfer from a 160 n.ml orbit, what Is the 
Earth-fixed longitude of deployment? (This Is the longitude of radar 
placement If the satellite is to be tracked from the ground during 
deployment.) 

ANS: 163° 51' 31” E 

D.11. SUMMARY FOR EARTH-SUN CALCULATIONS 
D.11.1. Symbols 

LMT = Local Mean Time 
GMT = Greenwich Mean Time 
LST = Local Sidereal T Ime 
GST = Greenwich Sidereal Time 
GHA = Greenwich Hour Angle 
RA = Right Ascension 
g = argument of perigee 
M = mean anomaly 
e = eccentricity 
£ = obliquity of plane of motion 

© = true anomaly 
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u = argument of latitude 

I s declination of object 

d = number of days since epoch 
X s earth fixed longitude, positive east 
^ E = earth fixed longitude east 

Xw = earth fixed longitude west 

LAT = earth fixed latitude 

p = central angle 

t = time 

to = time of epoch 

JD = Julian Day number 


D.11.2. Relations 


LMT 

= GMT+ XE 

LMT 

= GMT - Xw 

LST 

= QST+ Xe 

1ST 

= GST - X w 

QHA 

= GST 

u 

= g + 0 

1 

= sin-* (sin E sin u) 

RA(8un) 

= tan-' /X) 

X 

X 

= cos u 

y 

= cos e sin u 


Solar/Sidereal Conversion 

solar time Interval 
sidereal time Interval 


1.0027379093 


sidereal time interval 
solar time interval 


= 0.9972695664 


Arc/Tlm'3 Conversion ■ earth rotation relative to stars 


1 hour mean solar time = 15.014069 degrees of arc 
ArcTTIme Conversion - earth rotation relative to sun 

1 hour mean solar time = 15,00 degrees of arc 
Longitude of sunrise/sunset on earth 

-GHA(sun) ± cos-*(^°® ^ : 
sunset cos l,y„ cos LAT 

QHA (sun) = QHA (T) - RA (sun) 
Sunrise/sunset for spacecraft (s) 


RA (s) 


sunrise 

sunset 


RA (sun) ± cos-' ( 


Cos ^ - sin I 3 sin ^ 
cos i, cos i.,„ 


Earth Fixed Longitude for subsatellite point 


Xe = RA(s)-GHA(T) 


0.12. EVALUATION 
D.12.1. Learning Objectives 

The principal objectives of Section D were threefold; first, to introduce the 
measurement of relative position and time; second, to increase your problem 
solving capability in the conceptual range; and thirdly, to teach the methods for 
evaluating time related problems, such as launch window determination 

Having completed this section and related references, you should have an 
understanding of: 

• the relationship between various time measurements 

• the use of reference points, lines, and angles in the location 
of objects In space 

• the use of reference tables 

• the importance of earth-sun relations 
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D.12J. Outcome Meaeursit 


You should now be able to; 

1 . Convert time from one measurement set to another by use 
of formua or table. 

2. Convert time measurements to angular measure. 

3. Find the right ascensions, declinations, or other parameters 
of the sun and earth arbiters at specified times given appro- 
priate information. 

4. Find spacecraft longitude and longitudes of sunrise and 
sunset for earth and spacecraft. 

5. Calculate deploy windows. 

6. Consider advancing to changes of orbital plane and pertur- 
bations or going through Section D again, or asking lots of 
questions. 



APPENDIX A 


Selected data from 
JSCU262 

SPACE SHUTTLE ASTRODYNAMICAL CONSTANTS 

By B. F. Cockrell and Bruce Williamson 
Mathematical Physics Branch 


1.0 SUMMARY 

This document provides basic space Shuttle astrodynamic constants for use in 
mission planning and construction of ground and onboard software input loads. The data 
included here are provided to facilitate the use of consistent numerical values throughout 
the project. The document supersedes reference 1. 


2.0 INTRODUCTION 

The astrodynamic constants presented in this document are taken primarily from 
reference 2, "Natural Environmental and Physical Standards for Apollo and AAP". Those 
data have been expanded, and other data have been included when there was an indica- 
tion that additional information would be useful. The values of the constants in reference 
3 have been informally adopted by the Interplanetary Trajectory Committee. 


3.0 ASTRODYNAMIC CONSTANTS AND PARAMETERS 

3.1 Introduction 

This section presents the values and associated uncertainties of the constants and 
models used in trajectory prediction. A list of conversion factors and a description of the 
gravitational potential equation are also Included. 

Uncertainties (1 a ) are presented if available. For consistent conversion between units, 
more decimal digits are given for some quantities than are justified by the uncertainties. 
The values presented are in agreement with those adopted by NASA Headquarters (ref. 2). 


3.2 Constants And Parameters 

3.2.1. Angular Velocity of the Earth's Rotation 
With Respect to the Vernal Equinox 

The Earth’s angular rotational velocity with respect to a precessing equinox ( <^p) and 
Inertial equinox ( for OFT (calendar year 1979) based on references 8 and 16 is 
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uj p * ,7292115854918357 - 004 radians per second 
u), « ,7292115146459210 - 004 radians per second 
3,2,2. Speed of Sound Constant 

This is a derived constant used in the calculation of the speed of sound, The speed of 
sound is defined by: 


Cs 


' T 

y A ’m 
Mo 


Vi 

) 


where y = 1.4 and is the ratio of the specific heat of air at constant pressure to that at a 
constant volume (dimensimlers), 

Because , R*, and Mo are constants, they are combined for computational ease and the 
equation Is written; 



where K is the speed of sound constant 


K = 65,77035 ft/sec/ VlC 


3.2.3, Gravitational Potential Function 

The classical expression (or the gravitational potential V exerted at a point in space 
located at a distance r from the center of the attracting body of radius Rg and gravita- 
tional parameter p g Is given in equation as 

V Pnm (Sin 4 > ) c„„, cos (m X ) -f sin (m \) 
rn = 0 


<l> Is latitude (geodetic) 

X is longitude 

where c„^ and are the harmonic coefficients of the potential function and (sin 0 ) 
represents the associated Legendre functions of the first kind, of degree n and order m, 
Because sin 0 = z/r = u, where u is a direction cosine, the associated Legendre func- 
tions may be expressed as; 


V(r, <6 , X ) 


s 

n = 0 


r 
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TEXAS SOUTHERN UNIVERSITY 

HOUSTON, TIXA8 77004 


oci>artmcnt or f>Hvaies 


'November 7, 1981 


Mr. Robert Brown, Technical Monitor 
Lyndon B. Johnson Space Center 
National Aeronautics and Space Administration 
Houston, Texas 77058 

RE: NASA Grant NAG 9-6 

Dear Mr. Brown: 

Please find enclosed the communication prepared for 
you dated August 31, 1981, my anticipated date for trans- 
mittal of the final grant report. I did not transmit it 
then and thought that my conversation with you and 
Mr. Scarlett concerning the matter would suffice for the 
delay. Having offered to send the report minus one por- 
tion of typed material, you both thought that it should 
be sent complete and that a delay was acceptable. But it 
was an error on my part not to have had the fiscal report 
sent, since it is normally sent separately and was probab- 
ly of greater immediate significance. 

The small delay stretched itself and we now believe 
the report too long past due. In the enclosed material, 
therefore, is one part in a more or less complete but un- 
finished form. If it is acceptable to and agreeable with 
you, finished copies will be mailed to you and to the 
STIF as soon as I receive it in appropriate form (November 
30 projected) . 

It has been a pleasure working with you, and 1 wish 
to extend from Texas Southern University our hopes for 
continued cooperative efforts. 
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TEXAS SOUTHERN UNIVERSITY 

HOUSTON, TIXAS 77004 


Okfantmcnt or ^HvaiCB 


August 31, 1981 


TO: Robert H. Brown, Technical Monitor 

National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 

FROM: H. E. Blackwell 

Texas Southern University 

RE: NASA Grant NAG 9-6 


National Aeronautics and Space Administration (NASA) grant NAG 9-6 was 
awaided to Texas Southern 'University for work related to Space Shuttle flight 
and operation. The initial phase of the work was projected for one year be- 
ginning January 1, 1980. The project was divided into two areas: (1) a re- 

view, criticism, testing and preparation of course materials for training of 
space flight design paraprof essionals ; and (2) analysis and development of 
the kinematics and approximate flexibility models of the Remote Manipulator 
System (RMS) of the Space Shuttle. 

Course development proceeded from initial materials provided by the 
Johnson Space Center (JSC) . This work utilized principally the teaching and 
course material development expertise of TSU faculty in the design and docu- 
mentation of a syste',1 for the training of individuals in the mechanics of 
space flight design. 

A program of research was planned to delve into the subtler aspects 
of remote manipulator systems, systems similar to those that will be requir- 
ed for future applications in space. The initial investigations used the 
configuration planned for the Shuttle Transportation System. This approach 
was used to allow Texas Southern to build upon the substantial work already 
accomplished on this system, thereby gaining a significant jump toward ad- 
dressing the problems of more sophisticated systems. A significent aspect 
of this area of work was the transfer of computer programs used to investi- 
gate remote manipulator systems from JSC to TSU. These programs permit the 
rapid display and study of the control and response characteristics of Sys- 
tem manipulator and represent many years of analysis, coding, and testing. 
NASA was to have provided engineering data, algorithms used in remote con- 
trol, software documentation, and other information to aid TSU in building 
a remote manipulator capability. 

The first phase or initial grant period was listed as January 1 
through December 31, 1980, although the award was made several weeks after 
that beginning date. Funds for equipment (computer) came later in the year 
by way of a grant supplement. As a result of this and the incapacitation 
of the project director during the latter part of the year, a no-cost ex- 
tension of the grant period through June 1981 was obtained. Although sig- 
nificantly new developments in the RMS were not expected as a product of 
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this initial phase, the extension allowed completion of computer program trans- 
fer from JSC, along with other tasks, and was thought necessary to provide a 
more complete preparation for the second phase of this work. In May 19Sl, we 
were informed that funding and work force in this area at JSC had been reduced, 
and that funds for continued research at TSU along the lines of the initial 
grant could not be supported. 

A large portion of the initial work was designed to provide project fac- 
ulty with appropriate background, primarily the configuration and basic opera- 
tion of the Shuttle RMS along with the tools and skills necessary for a de- 
scription of that system. The RMS is a rather complex system having at that 
time a number of engineering problems remaining after years of work by a number 
of engineering groups. Our initial work consisted of learning tasks, movirig 
toward more specific problems of the RMS system and more creative approaches 
in developing such general systems which were to be accomplished in the second 
phase. 


Because most of the RMS work consisted of investigative tasks that were 
not new but oriented toward continuance, we will not detail that work. We 
are including, however, a report describing the RMS computer simulation pro- 
gram developed at JSC which we transferred to and worked with in our system: 
this includes work by E. L. Copeland of the Lockheed Engineering and Manage- 
ment Company whose assistance we acknowledge. One new approach on which we 
worked dealt with errors, in which a general analysis of error generation in 
r,i.e complete operation of the remote arm system was projected. Another con- 
sidered the spatl: 1 transformations necessary to describe end effector posi- 
tion and orientation: this approach used the Alpine concept in a move toward 

better control of time-response characteristics of the simulated system. Our 
principal reference for RMS study was documentation by SPAR Aerospace Products 
of Toronto, Canada. Study of Vought's Teleoperator Maneuvering System was an- 
ticipated as general development proceeded. 

Addressing the first area of the grant, course materials for training 
of space flight design paraprof essionals , we have included material labeled 
Introduction to Flight Planning (I). 
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The constants are for real time operations are; 

C,o = - 1082.7 X 10 * 

C 30 = 2.56 X 10 -* 

C, 3 = 1.58 X 10-' 

C 32 = 1.57 X 10-' 

Sjj = -.897 X 10-' 

all others are zero. 

(Higher order models are sometimes used for analysis and postmission trajectory 
reconstruction.) 


3.2. 3.1. Equatoriai Earth Radiuu (Gravitational) 

1 E.r. = .6378160000000000 + 007 ± .50000 + 001 m (ref. 2). 

1 E.r. = .2092572178477690 + 008 ± .16404 + 002 int. ft. 

1 E.r. = .3443930885519158 + 004 ± .26998-002 n. mi. 

1 E.r. = .9999992160754699 + 0C0 ± .78392-006 E.r. (MCC). 

3.2.3.2. Gravitational Parameter (GM^ = U ^ = M Earth) 

M e = .39860.2000000000 + 015 ± .40000 + 009 m^/secs (ref. 2). 

M e = .1407646853278542 + 017 ± .14126 + 001 (int. ft)3/sec2. 

M e = .6275027808522208 + 005 ± .62970 - 001 (n. mi.)3/sec2 

M e = .1990931661816326 + 002 ± .19979 - 004 (E.r.)s/hr2 (MCC). 

3.2.3.3. Mass of the Earth 


M 3 = .5973343323842350 + 025 KG 
M 3 = .1316896781981220 + 026 LBM 

3.2.4, Lunar Constants 


3.2.4. 1 . Earth-Moon Mass Ratio (ref. 2) 
M 3 /M 3 , = 81.3010 (±0.0010) 

3.2.4.2. Mean lunar Radius 


R 3 , = .1738090000000000 + 007 ± .70000 + 

R„, = .5702395013123360 + 007 ± .22966 + 

R„, = .9384935205183585 + 003 ± .37797 - 

R„ = .2725062772756741 + 000 ± .10975 - 


002 m (ref. 2 ). 

003 int. ft. 

001 n. mi. 

004 E.r. (MCC). 


3.2.4.3. Principal Axes 


a = .1738570000000000 + 007 ± .70000 + 002 m (ref. 2). 
a = .5703969816272966 + 007 ± .22966 + 003 int. ft. 
a = .9387526997840173 + 003 ± .37797 - 001 n. mi. 
a = .2725815340305558 + 000 ± .10975 - 004 E.r. (MCC). 
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b = .1738210000000000 + 007 ± .70000 -f 002 m (ref. 2). 

b = .5702788713910761 + 007 ± .22966 + 003 Int. ft. 

b = .9385583153347732 003 ± .37797 - 001 n. mi. 

b = .2725250914643945 + 000 ± .10975 - 004 E.r. (MCC). 

c = .1737490000000000 007 ± .70000 + 002 m (ref. 2). 

c = .5700426509186352 + 007 ± .22966 + 003 int. ft. 

c = .9381695464362851 + 003 ± ,37797 - 001 n, mi. 

c = .2724122083320720 + 000 ± ,10975 - 004 E.r. (MCC). 

where a is directed toward the center of the Earthj c is coincident with the Moon’s rota- 
tional aKis, and b Is perpendicular to a and c, 

3.2.4.4. Gravitational parameter for the Moon (GM,„ = m Moon) 

p Moon = .4902780000000000 -f 013 ± .60000 008 mVsec^ (ref. 2). 

(j Moon = ,1731400417087798 4- 015 4: ,21189 -4 010 (Int. ft)Vsec2. 

^ Moon = .7718260968373028 003 ± .94456 - 002 (n. mi.)3/sec!, 

p Moon = .2448838571715250 + 000 ± .29969 - 005 E.r.Vhr^ (MCC). 

3.2.5. General Constants 

3.2. 5.1. Astronomical Unit 

AU = .1495978930000000 4- 012 ± .50000 4 004 m (ref, 2), 

AU = .4908067355643045 4- 012 ± .16404 4- 005 Int. ft. 

AU = .8077640010799136 4 008 ± .26998 4 001 n, mi. 

AU .2345469159233102 4 005 ± ,78I'l93 - 003 E.r, (MCC), 

3. 2.5. 2. Velocity of Light In a Vacuum 

c = .2997925000000000 4 009 ± .30000 4 003 m/sec (ref. 2). 

c = ,9835711942257218 4 009 ± .98425 4 003 Int, ft/sec. 

c = .1618750000000000 4 006 ± .16199 4 000 n. ml./sec. 

c = ,1692105801590269 4 006 ± .16933 4 000 E.r./hr (MCC). 

3. 2.5. 3. Gravitational Parameters for the Sun 

M Sun = .1327124990000000 4 021 ± .15000 4 014 mVsec* (ref. 2). 

H Sun = .4686697671960388 4 022 ± .52972 4 015 (int. H)Vsec‘. 

H Sun = ,2089242635906454 4 011 ± .23614 4 004 (n. ml.)3/sec2. 

H Sun = .6628718533157138 4 007 ± ,74922 4 000 E.r.Vhr^ (MCC). 

3.3. EPHEMERiS TAPE SYSTEMS 

3.3.1. DE19 Tape 

The ephemerls tape system to be used for all missions Is provided by JPL and Is called 
the JPL Development Ephemerls Number 19 (DE19). For additional information, see 
references 7 and 8. 
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3,3,2, Tape System Conversion Factors 

The following values are to be used to convert DE19 units to Kilometers (refs, 7,8, and 

11 ). 

AU = 149 597 1)93 Km (scale factor for planetary ephemerldes) 

= 6378.1492 Km (scale factor for lunar ephemerldes) 

M ■* = 81,301 (ratio of Earth mass to Moon mass) 

GE = 398 601.2 kmVsec^ (gravitational parameter of the Earth) 


3.4. FISCHER EARTH MODEL 

The following constants describe the Fischer Earth model (1960), which Is used for the 
location of radar stations and other Earth surface features (ref, 2) 

3.4.1. Equatorial Earth Radius 

a = ,6378166000000000 + 007 m. 
a =* ,2092574146981627 + 008 Int, ft, 
a = .3443934125269978 + 004 n, ml. 
a = ,1000000156784906 + 001 E,r. (MCC), 

3.4.2. Flattening 

1 = flattening = 1 - b/a 
f = 1/298.30 = 0.3352329869259135 x 10 » 

3.4.3. Polar Earth Radius 

b = .6356784283607107 + 007 m, 
b = .2085559148165061 + 008 Int. ft. 
b .= .3432388922034075 + 004 n. ml. 
b = .9966478263900521 + 000 E.r. (MCC), 

3.4.4. Eccentricity of Ellipsoid 



e = .8181333401693114 - 001 
e» = 2f - f2 

= .6693421622965943 - 002 


4.0 EQUIVALENTS AND CONVERSION FACTORS 

1 Int, ft = .3048000000000000 + 000 m (exact) (ref, 2) 

1 n, ml. = ,1852000000000000 + 001 km (exact) (ref. 2) 

1 E.r. = 6378165000000000 = 004 km (exact for scaling) for MCC Internal use 
1 Ibm = .4535923700000000 + 000 kg (exact) (ref, 2) 



1 hr = 3600.000000000000 000 sec (exact) 

1 rad s: 180/n deg 
1 deg = 3600 arc sec 
1 km = .5399568034657235 + 000 n, ml. 

1 m = .3280839895013123 4- 001 Int. ft 
1 n. ml. = .60761 1&I85564304 4- 004 Int. ft. 

1 rad = .572957795t 308233 4- 002 deg 
1 deg = ,1745329251994329 - 001 rad 
1 kg = .2204622621848776 4- 001 Ibm 
1 Int. stat, ml. = 5280 ft (exact) 

1 Ibf = 32.174048556 (Int. ft/sec^) Ibm 
V = ,3141592653589793 4- 001 (ref, 2) 
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